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gravity 

gr 
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1 
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max 

maximum 
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1  Introduction 


Helically— grooved  heat  pipes  (HGHPs)  have  potential  applications  in  the  thermal 
management  of  rotating  equipment  such  as  aircraft  alternators,  large-scale  industrial 
electric  motors,  and  spinning  satellites.  In  two  recent  studies  (Klasing  et  al.,  1999; 
Thomas  et  al.,  1998),  the  performance  of  revolving  HGHPs  was  investigated.  It  was 
found  that  the  capillary  limit  increased  with  the  strength  of  the  acceleration  field 
perpendicular  to  the  heat  pipe  axis.  In  order  to  move  HGHPs  closer  to  application, 
knowledge  must  be  gained  concerning  the  sensitivity  of  the  capillary  limit  to  working 
fluid  fill  amount,  since  variations  in  the  fill  amount  are  inevitable  during  the  man¬ 
ufacture  of  these  devices.  Very  few  studies  were  available  concerning  the  effect  of 
working  fluid  fill  on  the  performance  of  axially-grooved  heat  pipes,  but  those  found 
have  been  outlined  below.  In  addition,  synopses  of  the  two  aforementioned  studies 
on  revolving  HGHPs  have  also  been  provided. 

Brennan  et  al.  (1977)  developed  a  mathematical  model  to  determine  the  per¬ 
formance  of  an  axially-grooved  heat  pipe  which  accounts  for  liquid  recession,  liquid- 
vapor  shear  interaction  and  puddle  flow  in  a  1-g  acceleration  environment.  The  model 
considered  three  distinct  flow  zones:  the  grooves  unaffected  by  the  puddle,  the  grooves 
that  emerge  from  the  puddle,  and  the  grooves  that  are  submerged  by  the  puddle.  The 
model  for  the  puddle  consisted  of  satisfying  the  equation  of  motion  for  the  puddle 
and  the  continuity  equation  at  the  puddle-groove  interface,  and  was  solved  by  a 
fourth-order  Runge-Kutta  integration  method  with  self-adjusting  step  sizes.  The 
assumptions  made  by  the  model  for  the  puddle  were  uniform  heat  addition  and  re¬ 
moval  with  a  single  evaporator  and  a  single  condenser  section,  and  one-dimensional 
Iflminflr  flow  in  the  puddle.  The  transport  capability  of  the  grooves  unaffected  by 
the  puddle  and  the  grooves  extending  beyond  the  puddle  were  approximated  by  a 
closed-form  solution  with  laminar  liquid  and  vapor  flow.  The  working  fluids  used  for 
the  experiment  were  methane,  ethane  and  ammonia.  Brennan  et  al.  (1977)  stated 
that  the  mathematical  model  agreed  well  with  the  experimental  data  for  ideally  filled 
and  overfilled  heat  pipes,  but  some  differences  were  noted  for  underfilled  heat  pipes. 
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In  general,  it  was  found  for  ideally  filled  heat  pipes  the  predicted  transported  heat 
was  higher  than  that  measured.  Also,  this  discrepancy  was  more  significant  for  lower 
operating  temperatures.  In  addition,  it  was  found  during  the  experiments  that  the 
maximum  transported  heat  increased  with  fill  volume. 

Vasiliev  et  al.  (1981)  performed  a  series  of  experiments  on  an  aluminum  axially- 
grooved  heat  pipe  which  was  overfilled  and  ideally  filled.  The  width  and  height  of 
the  grooves  were  w  =  0.123  mm  and  h  =  0.7  mm,  respectively,  with  an  overall  heat 
pipe  length  of  Lt  =  80.0  cm.  The  working  fiuids  were  acetone  and  ammonia.  Vasiliev 
et  al.  showed  that  the  temperature  difference  from  the  evaporator  to  the  adiabatic 
regions  increased  at  a  much  slower  rate  with  increasing  overfills.  This  was  attributed 
to  a  thin  film  of  liquid  emerging  from  the  overfill  pool  wetting  the  upper  grooves. 
Vasiliev  et  al.  stated  that  this  thin  film  was  lifted  over  the  grooves  by  capillary  forces 
due  to  microroughness  on  the  groove  surface.  A  mathematical  model  was  developed 
for  low  temperature  axially-grooved  heat  pipes  to  estimate  heat  pipe  performance 
for  0-g  and  1-g  applications.  The  mathematical  model  was  a  set  of  boundary-value 
problems  applied  to  each  groove  and  was  solved  by  a  numerical  iteration  method. 
The  model  was  based  on  pressure  balance  equations  and  mass  continuity  written  for 
a  single  groove.  The  temperature  of  the  vapor  in  the  adiabatic  region  was  an  input 
parameter,  and  the  vapor  pressure  gradient  was  assumed  to  be  one-dimensional.  In 
addition,  the  liquid-vapor  shear  stress  was  assumed  to  be  constant,  and  the  starting 
liquid  film  thickness  was  of  the  same  order  of  magnitude  as  the  groove  microroughness. 
Very  good  agreement  was  reported  between  the  mathematical  model  and  experimental 
transported  heat  results  for  ideally  filled  and  overfilled  heat  pipes  under  gravity. 

Thomas  et  al.  (1998)  presented  experimental  data  obtained  from  a  helically- 
grooved  copper  heat  pipe  which  was  tested  on  a  centrifuge  table.  The  heat  pipe 
was  bent  to  match  the  radius  of  curvature  of  the  table  so  that  uniform  transverse 
(perpendicular  to  the  axis  of  the  heat  pipe)  body  forces  field  could  be  applied  along 
the  entire  length  of  the  pipe.  The  steady-state  performance  of  the  curved  heat  pipe 
was  determined  by  varying  the  heat  input  (Qin  =  25  to  250  W)  and  centrifuge  table 
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velocity  (radial  acceleration  |ar|  =  0.01  to  10-g).  It  was  found  that  the  capillary  limit 
increased  by  a  factor  of  five  when  the  radial  acceleration  increased  from  |ar|  =  0.01  to 
6-g  due  to  the  geometry  of  the  helical  grooves.  A  model  was  developed  to  calculate 
the  capillary  limit  of  each  groove  in  terms  of  centrifuge  table  angular  velocity,  the 
geometry  of  the  heat  pipe  and  the  grooves,  and  the  temperature-dependent  working 
fluid  properties.  The  agreement  between  the  model  and  the  experimental  data  was 
satisfactory. 

Klasing  et  al.  (1999)  developed  a  mathematical  model  to  determine  the  operating 
limits  of  a  revolving  helically-grooved  straight  heat  pipe.  The  capillary  limit  calcu¬ 
lation  required  an  analysis  of  the  total  body  force  imposed  by  rotation  and  gravity 
on  the  liquid  along  the  length  of  the  helical  grooves.  The  boiling  and  entrainment 
limits  were  calculated  using  methods  described  by  Faghri  (1995).  It  was  found  that 
the  capillary  limit  increased  significantly  with  rotational  speed  due  to  the  helical  ge¬ 
ometry  of  the  heat  pipe  wick  structure.  The  maximum  heat  transport  was  found 
to  be  a  function  of  angular  velocity  and  tilt  angle  from  horizontal.  In  addition,  a 
TTiiniTrmm  value  of  angular  velocity  was  required  to  obtain  the  benefits  of  the  helical 
groove  geometry. 

The  first  objective  of  the  present  study  was  to  determine  the  sensitivity  of  the 
performance  of  revolving  HGHPs  to  the  working  fluid  fill  amount.  This  required  a 
precise  knowledge  of  the  geometry  of  the  heat  pipe  and  helical  grooves.  In  addition, 
a  precision  filling  station  was  constructed  and  calibrated  to  determine  the  uncertain¬ 
ties  involved  in  the  filling  procedure.  The  copper— ethanol  heat  pipe  was  tested  on 
a  centrifuge  table  at  Wright— Patterson  AFB  (AFRL/PRPG)  to  determine  the  cap¬ 
illary  limit,  thermal  resistance  and  evaporative  heat  transfer  coefficient  for  fill  ratios 
of  G  =  0.5,  1.0  and  1.5,  and  radial  accelerations  of  |ar|  =  0.01,  2.0,  4.0,  6.0,  8.0 
and  10.0-g.  The  second  objective  of  the  present  study  was  to  improve  the  existing 
analytical  capillary  limit  model  developed  by  Thomas  et  al.  (1998)  using  the  above- 
mentioned  geometric  measurements  and  by  using  improved  equations  for  the  working 
fluid  properties. 
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2  Determination  of  Heat  Pipe  Working  Fluid  Inventory 


The  objective  of  this  analysis  was  to  determine  the  working  fluid  inventory  of  a 
HGHP,  which  consists  of  the  mass  of  liquid  in  the  grooves  and  the  mass  of  vapor  in 
the  vapor  space.  Since  the  heat  pipe  is  a  closed  container  under  saturation  conditions, 
the  total  mass  of  working  fluid  in  the  heat  pipe  is  given  by 

Ks  GV„ 

=  7n„  +  m;  = - 1 - ^  (1) 

Vy  vi 

where  14s  aiid  Vgr  are  the  volumes  of  the  vapor  space  and  grooves,  Vy  and  vi  are  the 
speciflc  volumes  of  the  working  fluid  vapor  and  liquid,  and  G  =  Vi/Vgr  is  the  ratio  of 
the  volume  of  liquid  to  total  groove  volume.  The  volume  of  the  vapor  space  is  given 
by 


Ks  =  jDii,  +  v^(l  -  G)  (2) 

where  is  the  diameter  of  the  heat  pipe  vapor  space  and  Lt  is  the  total  heat  pipe 
length.  The  second  term  in  eqn.  (2)  accounts  for  the  increase  or  decrease  in  the  vapor 
space  volume  when  the  parameter  G  is  varied.  The  volume  of  the  grooves  is  given  by 


Vg[  —  Z/grTVgi-Agr 


(3) 


where  JVgr  is  the  total  number  of  grooves,  Ag^  is  the  groove  cross-sectional  area,  and 
Lgr  is  the  total  groove  length. 

A  cross-sectional  view  of  a  typical  helical  groove  in  the  experimental  test  article 
is  shown  in  Fig.  1.  The  cross-sectional  area  of  the  groove  is  given  by 


Agr  =  wh  +  -h^(tan  6i  +  tan  $2) 


(4) 


where  h  is  the  height  of  the  groove,  w  is  the  width  of  the  groove  base,  and  the  angles 
from  the  sides  of  the  groove  to  vertical  are  61  and  62.  The  total  length  of  each  groove 
is  given  by 


L 


gr 


—  Lt 


(5) 
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Figure  1:  Photomicrograph  of  the  helical  groove  geometry. 


The  radius  of  the  helix  is  given  by 

fh  —  2  (-^vs  +  h)  (6) 

The  helical  pitch  is  the  distance  through  which  the  helix  makes  one  revolution  around 
its  radius. 


_  27r(s  —  Si) 


(7) 


The  helix  angle  <f>  corresponds  to  s,  which  is  the  distance  traveled  along  the  centerline 


of  the  heat  pipe. 

In  order  to  calculate  the  working  fluid  inventory  for  the  HGHP,  measurements  of 
the  appropriate  geometric  parameters  were  made.  These  measurements  are  described 
in  the  following  section.  In  addition,  an  extensive  uncertainty  analysis  was  performed 
to  determine  the  uncertainties  of  both  the  measured  and  calculated  variables  used  in 


finding  the  working  fluid  inventory. 
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2.1  Geometric  Measurements 

The  physical  variables  given  in  eqn.  (4)  for  the  cross-sectional  area  of  the  grooves 
have  been  measured.  A  sample  of  the  HGHP  was  set  in  an  epoxy  resin  mold,  polished, 
and  examined  under  a  Nikon  (Model  Epiphot)  microscope  with  50  x  magnification. 
Buehler  Omnimet  computer  software  was  used  to  make  bitmap  pictures  of  ten  differ¬ 
ent  grooves  and  a  microscopic  calibration  scale.  These  pictures  were  then  analyzed 
using  Adobe  Illustrator  software  to  determine  the  geometric  values  shown  in  Fig.  1. 
Since  the  corners  at  the  top  of  the  land  between  grooves  were  not  well  defined,  a 
special  procedure  was  established  to  determine  the  geometry  of  the  grooves.  First, 
lines  were  drawn  along  the  bottom  and  sides  of  each  groove.  Then,  a  line  was  drawn 
across  the  bottom  of  the  land  between  grooves,  as  shown  in  Fig.  1.  This  line  was 
then  transposed  to  the  top  of  the  land.  The  intersection  between  this  line  and  the 
line  along  the  side  of  the  groove  was  defined  as  the  upper  corner  of  the  groove.  This 
procedure  was  also  performed  for  the  other  side  of  the  groove  to  establish  the  other 
corner.  Note  that  the  two  lines  along  the  land  tops  are  at  different  angles  due  to 
the  radius  of  curvature  of  the  heat  pipe  container.  The  angles  9i  and  62,  and  the 
height  and  width  of  the  groove  h  and  w  were  found  using  Adobe  Illustrator  software 
in  conjunction  with  the  microscopic  calibration  scale.  However,  the  resolution  of  the 
calibration  scale  was  20  /xm,  which  resulted  in  uncertainties  of  the  groove  area  on  the 
order  of  12%.  Therefore,  the  calibration  scale  was  evaluated  with  a  Nikon  (Model  R- 
14A)  optical  comparator  as  follows.  The  calibration  scale  was  mounted  in  the  optical 
comparator,  and  a  distance  of  300  fxm  indicated  by  the  calibration  scale  was  measured 
eight  times  as  shown  in  Table  1.  The  bitmap  picture  of  the  calibration  scale  was 
then  loaded  into  Adobe  Illustrator  and  measured  eight  times  using  the  Illustrator 
measurement  tool,  as  shown  in  Table  2  to  obtain  a  conversion  factor.  This  factor  was 
then  used  in  the  measurement  of  the  groove  height  and  width.  Table  3  shows  the 
raw  data  and  average  for  each  helical  groove  parameter  over  a  sample  of  ten  grooves. 

The  Nikon  optical  comparator  was  also  used  to  determine  the  vapor  space  diameter 
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Table  1:  Measurements  of  the  calibration  scale  by  the  optical  comparator. 


Trial 

Distance  (cm) 

1 

0.0301 

2 

0.0299 

3 

0.0299 

4 

0.0301 

5 

0.0300 

6 

0.0300 

7 

0.0301 

8 

0.0301 

Average 

0.030025 

Readability 

±0.0001 

Table  2:  Measurements  of  the  calibration  scale  by  the  Adobe  Illustrator  software. 


Trial 

Distance  (cm) 

1 

4.7272 

2 

4.7390 

3 

4.7272 

4 

4.7272 

5 

4.7390 

6 

4.7272 

7 

4.7272 

8 

4.7390 

Average 

4.731625 

Readability 

±0.0059 
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Table  3:  Statistical  results  for  the  helical  groove  geometry. 


Groove  Number 

h  (cm) 

w  (cm) 

ex 

02 

1 

0.03806 

0.03476 

15.08° 

14.08° 

2 

0.03838 

16.38° 

11.63° 

3 

0.03851 

0.03404 

15.89° 

13.88° 

4 

0.03958 

0.03462 

16.91° 

13.67° 

5 

0.03857 

0.03537 

15.69° 

14.77° 

6 

0.03791 

13.20° 

7 

0.03830 

0.03448 

15.48° 

13.74° 

8 

0.03777 

13.50° 

9 

0.03800 

14.82° 

10 

0.03806 

0.03432 

14.78° 

14.69° 

Average 

0.03831 

0.03445 

15.44° 

13.80° 

Standard  Deviation 

0.00052 

0.00080 

0.90° 

0.95° 

Table  4:  Statistical  results  for  vapor  space  diameter. 


Trial 

Dys  (cm) 

Trial 

Dys  (cm) 

1 

■Bll 

6 

1.3564 

2 

1.3638 

7 

1.3538 

3 

1.3640 

8 

1.3564 

4 

9 

1.3665 

5 

1.3564 

Average 

1.3586 

Standard  Deviation 

0.0049 

of  the  heat  pipe  container  sample.  The  cross  hairs  of  the  optical  comparator  were 
carefully  aligned  with  the  top  of  the  land  between  grooves  on  the  left  edge  of  the 
pipe.  The  comparator  table  was  then  moved  until  the  land  tops  on  the  right  edge 
of  the  pipe  were  aligned  with  the  cross  hairs.  The  diameter  of  the  heat  pipe  vapor 
space  was  the  distance  of  the  table  movement.  Table  4  shows  the  raw  data,  average 
and  standard  deviation  for  ten  vapor  space  diameter  measurements. 

The  helical  groove  pitch  was  found  using  a  vertical  milling  machine  and  an  angular 
displacement  transducer.  The  heat  pipe  container  material  was  originally  1  m  long. 
Approximately  one-half  was  used  to  form  the  heat  pipe,  and  the  other  half  was  used 


to  determine  the  pitch.  The  rotation  angle  (0  —  ^i)  and  the  corresponding  distance 
along  the  centerline  of  the  heat  pipe  (s  —  Si)  has  been  found  as  shown  in  Fig.  2(a). 

A  heat  pipe  holding  device  was  constructed  from  two  angle  aluminum  uprights 
mounted  to  the  table  of  a  Positool  vertical  milling  machine  (Model  S-1).  Precision 
alignment  blocks  were  attached  to  the  undersides  of  the  uprights  to  engage  one  of  the 
grooves  in  the  milling  machine  table  for  improved  alignment.  Nylon  bushings  were 
placed  in  the  uprights  to  center  both  the  heat  pipe  container  and  the  shaft,  which 
was  concentric  with  the  heat  pipe  container.  A  small  pin  was  made  from  a  1.58  mm 
(0.0625  in)  dowel  pin,  where  one  end  was  ground  to  0.26  mm  to  fit  in  the  base  of  the 
helical  groove.  This  sprung  pin  was  set  in  a  hole  in  the  shaft  where  it  engaged  one 
of  the  grooves,  as  shown  in  Fig.  2(6).  A  Trans-Tek  angular  displacement  transducer 
(Model  0605)  was  mounted  onto  another  piece  of  angle  aluminum.  A  vertical  6.35 
mm  (0.25  in)  dowel  pin  was  placed  in  the  angle  aluminum  to  align  with  the  angular 
displacement  transducer  shaft.  The  dowel  pin  was  held  by  a  collet  installed  in  the 
milling  machine  spindle  in  order  to  fix  the  location  of  the  displacement  transducer. 
The  shaft  of  the  transducer  was  linked  to  the  shaft  within  the  heat  pipe  by  three  set 
screws.  As  the  milling  machine  table  moved  the  pipe  over  the  stationary  shaft,  the 
pin  followed  the  helical  groove,  causing  the  shaft  to  rotate.  The  angular  displacement 
transducer  measured  this  rotation.  A  Hewlett-Packard  multimeter  (Model  3478A) 
was  used  to  measure  the  output  voltage  of  the  angular  displacement  transducer,  as 
shown  in  Fig.  3.  The  distance  of  the  table  movement  was  (s  -  si),  which  was  read 
from  the  milling  machine  display  unit.  The  transducer  output  voltage  was  measured 
over  10  cm  lengths  for  ten  different  groups  as  shown  in  Table  5.  Backlash  errors 
were  avoided  by  not  reversing  the  table  movement  while  taking  data.  The  pitch  was 
calculated  using  eqn.  (7)  at  a  point  in  the  center  of  each  10  cm  length.  The  helical 
pitch  versus  distance  along  the  heat  pipe  container  is  presented  in  Fig.  4. 

The  total  active  heat  pipe  length  was  measured  with  a  machinist’s  scale  and  a 
vernier  caliper.  The  active  length  was  the  total  container  length  minus  the  lengths  of 
the  two  end  caps  which  sealed  the  container. 
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1 .  Angular  Transducer  4.  Shaft 

2.  Milling  Machine  Spindle  5.  Heat  Pipe  Container 

3.  Milling  Machine  CoUet  6.  Milling  Machine  Table 


Figure  2:  Schematic  of  the  helical  pitch  measurement  technique:  (a)  Major  compo¬ 
nents;  (b)  Cross-sectional  view  of  sprung  pin  engaging  a  helical  groove. 
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Figure  3:  Photograph  of  the  pitch  measurement  setup. 


s  (cm) 

Figure  4:  Helical  pitch  measurement  data:  p  =  135.9  cm 
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Table  5:  Measurements  of  the  angular  displacement  voltage  versus  vertical  milling 
machine  travel. 


Table 
Movement 
s  (cm) 

Transducer 

Voltage 

Table 
Movement 
s  (cm) 

Transducer 

Voltage 

Table 
Movement 
s  (cm) 

Transducer 

Voltage 

0.0 

0.9246 

1.0 

0.9992 

2.0 

1.1477 

10.0 

2.3569 

11.12 

2.5243 

12.0 

2.6589 

20.0 

3.8465 

21.0 

3.9982 

22.0 

4.1492 

30.0 

5.3273 

31.0 

5.478 

32.0 

5.6296 

40.0 

6.776 

41.0 

6.898 

42.0 

7.0338 

40.0 

6.7562 

41.0 

6.8962 

42.0 

7.0164 

30.0 

5.363 

31.0 

5.5093 

31.998 

5.6641 

20.0 

3.8876 

21.0 

4.0316 

22.0 

4.1789 

10.0 

2.38455 

10.999 

2.5328 

12.0 

2.684 

0.0 

0.86335 

1.0 

0.9996 

1.999 

1.1506 

3.001 

1.3068 

4.0 

1.4679 

5.001 

1.6183 

13.0 

2.792 

14.0 

2.9442 

15.001 

3.1004 

23.0 

4.2874 

24.0 

4.4335 

25.006 

4.5868 

33.0 

5.7794 

34.0 

5.9242 

35.0 

6.0517 

43.0 

7.1788 

33.999 

5.9354 

35.0 

6.0368 

43.0 

7.2109 

24.0 

4.4586 

25.0 

4.6110 

33.0 

5.817 

14.0 

2.9745 

15.0 

3.1278 

22.998 

4.3179 

4.0 

1.471 

5.0 

1.6219 

13.0 

2.8270 

3.0 

1.3147 

1.7673 

1.9100 

8.0 

2.0462 

3.2564 

17.0 

3.4052 

18.0 

3.5462 

4.7381 

4.8980 

28.0 

5.0431 

6.1890 

37.0 

6.3369 

38.0 

6.4862 

6.2184 

37.0 

6.3304 

38.0 

6.4825 

27.0 

4.9256 

28.0 

5.0670 

15.999 

17.0 

3.4299 

18.0 

3.5739 

6.0 

1.7740 

1.9176 

8.0 

2.0638 

2.2038 

3.6977 

5.1804 

6.6277 

3.7255 

2.2243 
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2.2  Uncertainty  Analysis 


Using  the  analysis  given  by  Miller  (1989),  the  root-sum-square  uncertainties  for 
the  groove  cross-sectional  area,  helical  pitch,  helix  radius,  groove  length,  groove  vol¬ 
ume,  vapor  space  volume,  and  total  mass  of  the  working  fluid  inventory  are  as  follows. 
Groove  cross-sectional  area: 


Helical  pitch:  The  rotation  angle  of  the  angular  displacement  transducer  (in  arc 
degrees)  was  given  by  the  following  calibration  equation 

(0  -  <^i)  =  17.986(V  -  Vi)  (9) 


where  V  is  the  output  voltage.  Using  this  expression  results  in  the  following  equation 


for  helical  pitch. 


20.015(s  -  si) 


(V-V:) 

Therefore,  the  uncertainty  of  the  helical  pitch  measurement  is 

f  / 20.015  A(s- Si)  V  ,  f20.015(s-si)  A(V-Vi) 
(V-Vi)  )^[  (V-Vi)2 
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Groove  length: 


Groove  volume: 


Vapor  space  volume: 


Total  mass  of  working  fluid  inventory: 


In  general,  the  uncertainties  of  the  measured  values  were  comprised  of  the  read¬ 
ability  and/or  uncertainty  of  the  measurement  tool  and  the  standard  deviation  found 
from  several  measurements.  To  be  conservative,  the  greater  of  the  readability  or  the 
uncertainty  of  the  measurement  tool  was  used  in  the  calculations.  The  values  of  stan¬ 
dard  deviation  for  h,  w,  62  and  Dys  are  given  in  Tables  3  and  4.  The  readability 
or  uncertainty  of  each  instrument  used  in  the  measurements  is  provided  in  Table  6, 
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Table  6:  Readability/uncertainty  values  for  the  instruments  used  in  determining  the 
working  fluid  inventory. 


Instrument 

Readability/ 

Uncertainty 

Variables 

Adobe  Illustrator  Angle 
Tool 

±0.01° 

$1,  02 

Adobe  Illustrator  Dis¬ 
tance  Tool 

±0.000037  cm 

h,  w 

Optical  Comparator 

±  0.0001  cm 

D-vs 

Microscopic  Calibration 
Scale 

±  0.0001  cm 

h,  w 

Multimeter 

±  0.00455  V 

V-Vi 

Angular  Displacement 
Transducer 

±0.0083  V 

V-Vi 

Vertical  Milling  Machine 

±  0.001  cm 

S  —  $i 

Machinist’s  Scale 

±0.079  cm 

Lt 

Vernier  Caliper 

±0.0025  cm 

Lt 

which  also  indicates  the  variables  found  by  each  instrument.  The  groove  height  and 
width  uncertainties  included  the  readabilities  of  the  optical  comparator,  the  micro¬ 
scopic  calibration  scale  and  the  Adobe  Illustrator  distance  tool  (Tables  1  and  2). 
The  measured  and  calculated  uncertainties  for  all  geometric  variables  presented  are 
shown  in  Table  7. 

2.3  Working  Fluid  Properties 

A  literature  survey  was  completed  to  determine  the  speciflc  volumes  of  ethanol 
vapor  and  liquid  at  various  saturation  temperatures,  as  shown  in  Fig.  5.  This  infor¬ 
mation  was  needed  to  determine  the  total  mass  and  uncertainty  of  the  working  fluid 
inventory  mt±Amt  [eqns.  (1)  and  (16)].  While  existing  texts  report  these  properties 
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Figure  5:  Specific  volume  of  ethanol  versus  temperature:  (a)  Saturated  liquid;  (6) 
Saturated  vapor. 
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Table  7:  The  geometric  variable  values  associated  with  the  working  fluid  inventory. 


Measured  Values 

h 

0.03831  ±  0.00076  cm 

w 

0.03445  ±  0.0010  cm 

01 

15.44®  ±  0.91® 

62 

13.80®  ±  0.96® 

V- Vi 

1.474  ±  0.064  Volt 

S  —  Si 

10.000  ±  0.023  cm 

Dys 

1.359  ±  0.005  cm 

Lt 

43.8  ±  0.084  cm 

Calculated  Values 

Agr 

1.703  X  10“®  ±  6.0  X  10"^  cm^ 

V 

135.8  ±  5.9  cm 

Th 

0.6992  ±  0.0025  cm 

Lgr 

43.82  ±  0.84  cm 

3.73  ±  0.13  cm^  ' 

(Faghri,  1995;  Peterson,  1994,  Lide  and  Kehiaian,  1994;  Carey,  1992;  Schlunder,  1983; 
Ivanovskii  et  al.,  1982),  it  was  found  that  most  simply  referred  to  previous  sources. 
Therefore,  the  data  shown  in  Fig.  5  represent  information  gathered  from  primary 
sources  that  cannot  readily  be  traced  further.  In  Fig.  5(a),  the  available  data  for 
the  specific  volume  of  liquid  in  the  range  of  Tgat  =  0  to  100°C  are  relatively  scat¬ 
tered.  Vargaftik  (1975)  stated  that  the  ethanol  used  was  96%  pure  by  volume,  with 
water  making  up  most  of  the  other  4%.  Ethanol  is  agressively  hygroscopic,  so  special 
procedures  are  required  for  further  purification  as  outlined  by  Timmermans  (1950) 
concerning  anhydrous  ethanol.  Since  the  data  by  Timmermans  (1950)  and  TRC 
(1983)  are  nearly  coincident,  it  is  believed  that  the  data  reported  by  TRC  (1983)  are 
also  for  anhydrous  ethanol.  Dunn  and  Reay  (1978)  do  not  provide  information  con¬ 
cerning  purity.  Therefore,  the  Vargaftik  (1975)  data  and  the  Dunn  and  Reay  (1978) 
data  have  been  discarded  in  Fig.  5(a).  In  Fig.  5(6),  the  deviation  of  the  Dunn  and 
Reay  (1978)  data  for  the  specific  volume  of  vapor  is  significant.  Therefore,  the  Dunn 
and  Reay  (1978)  data  has  been  discarded  in  Fig.  5(6).  Polynomial  curve  fits  from 
0  <  Tsat  <  100°C  have  been  obtained  for  the  data  shown  in  Figs.  5(a)  and  5(6)  for 
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Table  8:  The  calculated  total  mass  of  the  working  fluid  inventory. 


G 

int  (g) 

Amt  (g) 
Calculated 
[eqn.  (16)] 

Amd  (g) 
Filling 
Station 

0.5 

1.47 

±  3.6% 

±  5.0% 

1.0 

2.92 

±  3.7% 

±  2.9% 

1.5 

4.38 

±  3.6% 

±  1.9% 

the  specific  volumes  of  liquid  and  vapor  ethanol.  These  curve  fits  have  been  evaluated 
at  room  temperature  to  determine  the  proper  values  to  be  used  in  the  uncertainty 
analysis,  since  the  heat  pipe  was  filled  at  room  temperature.  Information  concerning 
the  uncertainty  of  the  original  data  was  not  available.  Therefore,  the  uncertainties  of 
these  properties  have  been  estimated  to  be  the  maximum  variance  of  the  data  from 
the  curve  fits  {Avi=  3.5  xlO“®  m^/kg,  Avy  =  0.39  m^/kg).  The  specific  volumes 
of  liquid  and  vapor  ethanol  (m®/kg)  as  functions  of  saturation  temperature  (®C)  are 
shown  below  for  the  range  0  <  Tgat  <  100°C 

Vi  =  exp(ao  +  CiTsat  +  0,2'lfa.t  +  Os^sat  +  ®4^sat)/1000  (17) 


Vy  —  exp(6o  +  5lTsat  +  627’sat  +  +  64Tsat)/1000  (18) 


where  the  coefficients  are 


oo  =  0.2153 
Cl  =  1.049  X  10-3 
02  =  -1.345  X  10-3 
as  =  2.025  X  lO-® 

04  =  -5.474  X  10-^1 


bo  =  10.35 

=  -6.375  X  10-3 

62  =  1.735  X  10-^ 

63  =  5.714  X  10-'^ 

64  =  -6.003  X  10-3 


The  total  mass  of  the  working  fluid  inventory  m*  and  the  associated  uncertainty  Amt 
for  the  range  of  fill  values  are  given  in  Table  8. 
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3  Heat  Pipe  Filling  Station 


3.1  Description 

A  filling  station  has  been  constructed  which  is  capable  of  placing  a  low-temperature 
working  fluid  (i.e.,  water,  ethanol,  methanol)  into  a  heat  pipe  without  also  introduc¬ 
ing  ambient  air.  A  schematic  of  this  station  is  shown  in  Fig.  6.  The  station  consisted 
of  a  manifold  of  valves  and  interconnecting  stainless  steel  tubing,  a  working  fluid 
reservoir,  a  dispensing  burette,  a  vacuum  pump,  and  a  container  of  compressed  dry 
nitrogen  gas.  Previous  experience  with  filling  stations  showed  that  long  runs  of  hor¬ 
izontal  tubing  could  cause  significant  filling  errors  due  to  vapor  bubbles  within  the 
tubing.  To  address  this  problem,  the  manifold  was  constructed  such  that  the  inter¬ 
connecting  tubing  runs  were  very  short  (on  the  order  of  2  cm).  In  addition,  the  tubes 
which  intersect  the  main  vertical  tube  between  valves  2  and  5  (Fig.  6)  were  offset 
from  each  other  and  ran  at  a  diagonal  from  the  main  tube.  Again,  the  purpose  of  this 
design  was  to  reduce  the  possibility  of  vapor  bubbles  adhering  to  the  tubing  walls, 
thus  causing  errors  in  the  fill  amount.  However,  it  is  likely  that  some  vapor  still  does 
adhere  to  the  tubing,  so  certain  procedures  were  carried  out  during  filling  to  eject  as 
much  vapor  as  possible.  For  instance,  the  1  psig  relief  valve  over  valve  1  was  cycled 
on  and  off  several  times.  In  addition,  valves  2  and  5  were  cycled  on  and  off  while 
noting  the  meniscus  displacement  within  the  dispensing  burette.  If  the  meniscus  was 
displaced  more  than  0.06  cm^,  vapor  was  probably  trapped  within  the  valve.  The 
valve  in  question  was  then  cycled  until  the  bubble  was  ejected. 

To  fill  the  heat  pipe,  the  container  was  first  evacuated  to  a  pressure  of  10“®  Torr 
using  a  Varian  turbomolecular  vacuum  pump  (Model  V450).  The  sealed  pipe  was  then 
connected  to  the  filling  station  at  valve  5.  The  working  fluid  was  frozen  and  thawed 
repeatedly  to  reduce  the  amount  of  dissolved  air  within  the  fluid.  The  entire  filling 
station  was  then  evacuated  by  a  roughing  pump,  except  the  working  fluid  reservoir. 
After  evacuation,  the  liquid  working  fluid  was  drawn  up  into  the  dispensing  burette 
and  into  all  interconnecting  tubing.  After  noting  the  height  of  the  meniscus,  the 
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Working  Fluid 
Reservoir 


Heat  Pipe 


Figure  6:  Schematic  of  the  heat  pipe  filling  station. 


desired  amount  of  working  fluid  was  metered  into  the  heat  pipe  by  carefully  opening 
the  heat  pipe  fill  valve  8.  The  difference  in  height  of  the  liquid  column  was  related  to 
the  dispensed  mass  of  working  fluid.  A  full  description  of  the  step-by-step  process  is 
provided  in  Section  9.1. 

3.2  Uncertainty  Analysis  of  the  Dispensed  Msiss  of  Working  Fluid 

During  initial  testing  of  the  filling  station,  it  was  found  that  the  mass  of  working 
fluid  dispensed  into  the  heat  pipe  container  was  different  than  what  was  indicated 
by  the  dispensing  burette.  Therefore,  a  rigorous  calibration  of  the  filling  station  was 
undertaken  to  determine  a  correlation  between  the  change  in  volume  read  by  the 
dispensing  burette  and  the  change  in  mass  of  a  receiving  burette  attached  at  valve 
5,  which  was  measured  using  a  precision  scale  (Mettler  Model  PC4400).  Figure  7(a) 
shows  eight  sets  of  data  taken  at  seven  different  fill  levels  during  the  filling  station  cal¬ 
ibration  process  and  a  linear  curve  fit  to  the  data.  Figure  7(6)  indicates  the  maximum 
residual  of  the  data  from  the  linear  curve  fit.  As  expected,  the  residual  percentage 
decreased  as  the  dispensed  working  fluid  amount  increased.  The  maximum  residual 
was  added  to  the  readability  of  the  dispensing  burette  (±  0.02  cm®)  to  determine 
the  total  uncertainty  of  the  working  fluid  inventory  dispensed  by  the  heat  pipe  filling 
station  Am^  as  a  function  of  the  groove  percentage  (Table  8). 
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Figure  8:  Experimental  setup  for  testing  the  helically-grooved  heat  pipe. 


4  Experimental  Setup 

The  purpose  of  the  experiment  was  to  examine  the  steady-state  performance  of  a 
helically-grooved  copper-ethanol  heat  pipe  under  various  heat  inputs  and  transverse 
body  force  fields  using  a  centrifuge  table  located  at  the  Air  Force  Research  Laboratory 
located  at  Wright-Patterson  AFB  (AFRL/PRPG).  Specifically,  the  amount  of  work¬ 
ing  fiuid  was  varied  {G  =  0.5, 1.0  and  1.5)  to  determine  the  effects  of  under/overfilling 
on  the  capillary  limit,  thermal  resistance  and  evaporative  heat  transfer  coefficient  of 
the  HGHP.  To  ensure  uniform  radial  acceleration  fields  over  the  length  of  the  heat 
pipe,  the  pipe  was  bent  to  match  the  radius  of  curvature  of  the  centrifuge  table  {R  = 
1.22  m)  as  shown  in  Fig.  8.  Physical  information  concerning  the  heat  pipe  is  given 
in  Table  9.  It  should  be  noted  that  the  total  helix  angle  was  very  small:  Each 
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Table  9:  Helically-grooved  heat  pipe  specifications. 


Working  fiuid 

Ethanol 

Working  fiuid  charge 

mt  =  1.47,  2.92  and  4.38  g 

Evaporator  length 

Lg  =  15.2  ±  0.16  cm 

Adiabatic  length 

La  =■  8.2  ±  0.16  cm 

Condenser  length 

Lc  =  15.2  ±  0.16  cm 

Tube  outside  diameter 

Do  =  1-588  ±  0.005  cm 

Tube  wall  thickness 

ty,  =  0.0757  cm 

Radius  of  curvature 

R  =  1.22  m 

Wall/ wick  materials 

Copper 

Wick  structure 

Helical  grooves 

Number  of  Grooves 

Ngr  =  50 

Heater  element 

Nichrome  heater  tape 

Fill  valve 

Nupro  B-4HW  bellows  valve 

Calorimeter 

1/8  in.  OD  coiled  copper  tubing 

groove  rotated  through  an  angle  of  approximately  2.03  rad  (116  arc  degrees)  over  the 
length  of  the  pipe.  The  heat  pipe  was  mounted  to  a  platform  overhanging  the  edge 
of  the  horizontal  centrifuge  table.  This  allowed  the  heat  pipe  to  be  positioned  such 
that  the  radius  of  curvature  was  equivalent  to  the  outermost  radius  of  the  centrifuge 
table.  Insulative  mounting  blocks  were  used  to  ensure  that  the  heat  pipe  matched  the 
prescribed  radius  as  closely  as  possible.  The  horizontal  centrifuge  table  was  driven 
by  a  20-hp  dc  motor.  The  acceleration  field  near  the  heat  pipe  was  measured  by 
a  Columbia  Research  Laboratories  (Model  SA-307-TX)  triaxial  accelerometer.  The 
acceleration  field  at  the  centerline  of  the  heat  pipe  radius  was  calculated  from  these 
readings  using  a  coordinate  transformation. 

A  Clayborn  Labs  (Model  J-16-2)  pressure-sensitive  nichrome  heater  tape  with  an 
aluminized  backing  was  uniformly  wound  around  the  circumference  of  the  evaporator 
section  for  heat  input.  Power  was  supplied  to  the  heat  pipe  evaporator  section  by 
a  Kepco  (Model  ATE150-7m)  power  supply  through  power  slip  rings  to  the  table. 
While  the  current  reading  could  be  made  directly  using  a  precision  ammeter,  the 
voltage  across  the  electric  heater  had  to  be  measured  on  the  rotating  table  because 
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of  the  voltage  drop  between  the  control  room  and  the  table.  Therefore,  the  voltage 
at  the  heater  was  obtained  through  the  instrumentation  slip  ring  assembly  and  read 
by  a  Hewlett  Packard  (Model  3478A)  precision  multimeter. 

The  calorimeter  consisted  of  a  length  of  1/8  in.  OD  copper  tubing  wound  tightly 
around  the  condenser  section.  The  size  of  the  tubing  was  chosen  to  be  small  to 
minimize  the  eflfects  of  acceleration  on  the  performance  of  the  calorimeter.  Thermal 
grease  was  used  between  the  heat  pipe  and  the  calorimeter  to  decrease  contact  resis¬ 
tance.  Type  T  thermocouples  were  inserted  through  brass  T-branch  connectors  into 
the  coolant  inlet  and  exit  streams,  and  a  Signet  (Model  8511)  high-resolution  digital 
flow  meter  was  used  to  measure  the  mass  flow  rate  of  the  coolant  (50%  by  mass  ethy¬ 
lene  glycol/water  mixture).  The  mass  flow  rate  was  controlled  using  a  high-pressure 
booster  pump,  which  aided  the  low-pressure  pump  in  the  recirculating  chiller.  The 
temperature  of  the  coolant  was  maintained  at  a  constant  setting  by  the  recirculating 
chiller  (Tew  =  30®C).  Coolant  was  delivered  to  the  centrifuge  table  via  a  double-pass 
hydraulic  rotary  coupling.  The  mass  flow  rate  was  constant  for  all  experiments.  Val¬ 
ues  of  the  specific  heat  of  ethylene  glycol/water  mixtures  have  been  obtained  from 
ASHRAE  (1977),  which  were  in  terms  of  percent  ethylene  glycol  by  weight  and  tem¬ 
perature.  The  average  temperature  between  the  calorimeter  inlet  and  outlet  was  used 
to  evaluate  the  specific  heat.  The  specific  heat  did  not  vary  appreciably  during  testing 
since  it  is  a  weak  function  of  temperature. 

The  Signet  high-resolution  digital  flow  meter  was  calibrated  using  a  Mettler 
(Model  PC4400)  scale  and  a  stop  watch.  The  temperature  of  the  ethylene  glycol 
solution  was  set  to  Tin  =  30°C,  which  was  measured  at  the  calorimeter  inlet  location. 
The  temperature  was  monitored  by  Viewdac  data  acquisition  software.  The  flow  was 
then  allowed  to  fill  a  graduated  cylinder  for  60  s,  after  which  the  mass  of  the  ethylene 
glycol  solution  was  recorded.  This  was  repeated  several  times  at  diflPerent  flow  rates. 
The  signal  coming  from  the  flow  meter  was  displayed  by  a  Signet  8511  Compak  flow 
transmitter  with  the  K-factor  set  to  2629.  The  results  of  the  calibration  are  given  in 
Table  10.  The  time  and  mass  values  shown  in  the  table  were  averaged  from  the  raw 
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Table  10:  Averaged  flow  rates  at  different  Signet  meter  readings. 


Signet  meter 

reading 

0.14 

60.0 

185.64 

3.10  ±  0.05 

0.15 

60.0 

207.06 

3.45  ±  0.06 

0.16 

60.0 

211.80 

3.53  ±  0.06 

data.  The  percentage  of  ethylene  glycol  in  the  coolant  was  measured  in  the  following 
manner.  The  specific  gravity  of  the  coolant  was  found  by  mixing  530  g  of  water  to  530 
g  of  pure  ethylene  glycol,  which  was  measured  with  a  Cole-Parmer  (Model  08292-16) 
hydrometer.  The  mass  of  the  water  and  ethylene  glycol  each  was  measured  with  a 
Mettler  (Model  PC4400)  scale.  The  hydrometer  measurement  was  within  1.5%  of  the 
ratio  of  masses  obtained  from  weighing  the  water  and  ethylene  glycol  separately  be¬ 
fore  mixing  them  and  the  ASHRAE  (1977)  values.  The  percentage  of  ethylene  glycol 
was  measured  periodically  during  testing  to  ensure  that  the  mixture  did  not  change. 

Heat  pipe  temperatures  were  measured  by  Type  T  surface-mount  thermocouples, 
which  were  held  in  place  using  Kapton  tape.  Mounting  locations  for  the  thermocou¬ 
ples  are  shown  in  Fig.  9.  A  short  unheated  length  next  to  the  evaporator  end 
cap  was  instrumented  with  thermocouples  specifically  for  the  detection  of  dryout  in 
the  evaporator  section  and  accurate  thermal  resistance  measurements.  In  addition, 
groups  of  four  thermocouples  were  arranged  around  the  circumference  of  the  heat  pipe 
at  stations  in  the  evaporator  section  for  local  heat  transfer  coeflScient  information. 
Temperature  signals  were  conditioned  and  amplified  on  the  centrifuge  table.  These 
signals  were  transferred  off  the  table  through  the  instrumentation  slip  ring  assembly, 
which  was  completely  separate  from  the  power  slip  ring  assembly  to  reduce  electronic 
noise.  Conditioning  the  temperature  signals  prior  to  leaving  the  centrifuge  table 
eliminated  difficulties  associated  with  creating  additional  junctions  within  the  slip 
ring  assembly.  Temperature  and  acceleration  signals  were  collected  using  a  Keithley 


26 


Inboard 


T 

TC# 

jc(mm) 

TC# 

(mm) 

a  -  25.4  mm  > 

no 

22.2 

70 

218. 

Lg  =  152.4  mm  j 

\  2IT0B 

54.0 

81 

238. 

La  =  82.3  mm  i  i 

io  3IT0B 

92.1 

90 

259. 

Lf.  =  152.4  mm  \ 

/  4IT0B 

130. 

lOIO 

435. 

h  =  25,4  mm  \ 

/  5IT0B 

168. 

11 

Coolant  Inlet 

End  caps:  9.53  mm 

B 

61 

198. 

12 

Coolant  Outlet 

Figure  9:  Thermocouple  locations  and  relevant  lengths. 


Table  11:  Data  acquisition  mapping. 


Data 

Acquisition 

Channel 

Relay 

Number 

Switch 

ON 

Switch 

OFF 

TCO 

- 

l-I 

l-I 

TCI 

- 

1-0 

1-0 

TC2 

- 

lO-I 

lO-I 

TC3 

- 

11 

11 
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(Model  500A)  DAQ  conditioner  and  a  personal  computer  with  Viewdac  data  logging 
software.  Since  a  shortage  of  thermocouple  channels  existed  on  the  centrifuge  table, 
a  series  of  three  electrical  relays  were  engaged  to  read  one  set  of  thermocouples,  and 
disengaged  to  read  the  other  set,  as  shown  in  Table  11.  Also  shown  in  this  table  is 
the  mapping  between  the  thermocouple  reference  name  and  data  acquisition  channel 
number. 

Since  the  heat  pipe  assembly  was  subjected  to  air  velocities  due  to  the  rotation  of 
the  table  (up  to  11  ra/s  =  25  mi/hr),  efforts  were  made  to  reduce  convective  heat  losses 
from  the  exterior  of  the  heat  pipe.  A  thin-walled  aluminum  box  was  fabricated  to  fit 
around  the  heat  pipe.  Ceramic  wool  insulation  was  placed  inside  the  box  and  around 
the  heat  pipe  through  three  small  doors  on  the  top  of  the  box.  This  insulation/box 
arrangement  provided  an  effective  barrier  to  convective  losses  from  the  heat  pipe  to 
the  ambient. 

The  helically-grooved  copper-ethanol  heat  pipe  was  tested  in  the  following  manner. 
The  recirculating  chiller  was  turned  on  and  allowed  to  reach  the  setpoint  temperature, 
which  was  measured  at  the  calorimeter  inlet  (Tin  =  30°C).  The  coolant  flow  rate  was 
set  to  rhcw  =  3.45  g/s  (Signet  reading  of  0.15).  The  centrifuge  table  was  started 
from  the  remote  control  room  at  a  slow  constant  rotational  speed  to  prevent  damage 
to  the  power  and  instrumentation  slip  rings.  In  this  case,  the  radial  acceleration 
was  less  than  |5r|  <  0.01-g.  In  all  cases,  the  centrifuge  table  rotated  in  a  clockwise 
direction  as  seen  from  above.  Power  to  the  heater  was  applied  (Qjn  =  10  W)  and  the 
heat  pipe  was  allowed  to  reach  a  steady-state  condition,  which  was  determined  by 
monitoring  the  evaporator  temperature  detected  by  thermocouple  2-1  until  changes  of 
less  than  0.1°C  over  2  min.  were  noted.  The  power  to  the  heater  was  then  increased 
to  Qin  =  20  W  and  again  the  heat  pipe  was  allowed  to  reach  a  steady-state  condition. 
This  was  repeated  until  the  maximum  allowable  evaporator  temperature  was  reached 
(heater  .max  =  100°C)  or  the  maximum  power  input  was  reached  (<3in,max  =  200  W). 
After  all  data  had  been  recorded  the  power  to  the  heater  was  turned  off,  and  the  heat 
pipe  was  allowed  to  cool  before  shutting  down  the  centrifuge  table. 
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The  power  input  to  the  electric  heater  was  calculated  using  the  average  of  five 
amperage  and  voltage  readings  taken  during  each  test.  Also,  a  correction  was  made 
for  the  losses  generated  by  a  short  length  of  wire  leading  to  the  heater.  Therefore, 
the  heat  input  is 

Qin  =  IV-l\R^re)  (19) 

The  heat  removed  from  the  heat  pipe  by  the  calorimeter  is  given  by 

Qt  =  mCp (Tout -Tin)  (20) 


The  thermal  resistance  of  the  heat  pipe  was  defined  as  the  ratio  of  the  temperature 
drop  across  the  heat  pipe  to  the  transported  heat 

rj _ Teec  —  Tcec  /n-i  \ 

•^th  ^  '  ' 

Uout  “ -iinj 

where  T^ec  is  the  evaporator  end  cap  temperature  measured  by  TCl-I,  and  Tcec  is  the 
condenser  end  cap  temperature  measured  by  TCIO-L  The  heat  transfer  coefficient  in 
the  evaporator  section  is 


where  T^j  represents  the  external  wall  temperatures  measured  in  the  evaporator  sec¬ 
tion.  Using  the  analysis  given  by  Miller  (1989),  the  uncertainties  for  the  heat  input, 
transported  heat  thermal  resistance  and  heat  transfer  coefficient  are 

AQi.  =  JliV  -  2IR^„)  A/J'  +  (lAVf  +  {PAR^f  (23) 


A<5t  =  ^ (Tout  —  2in)  Aril]^  +  [lh  (Tout  —  Tin)  ACp]^  +  (mCpATin)  +  {ihCp^Tcnt) 

(24) 


A/2th  = 


r  (Teec  —  Tcec)  ACp  ^  (Teec  ~  Tgec)  ATjn  ^  (Teec  Tgec)  ATput 

(Tout- Tin)  J  [mCp  (Tout  -  Ti„)"J  mCp  (Tout  -  Tin)'  . 


1 

2 


(25) 
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Table  12:  Maximum  uncertainties  of  measured  and  calculated  values. 


Measured  Values 

Mass  flow  rate 

Am  =  ±  0.05  g/s 

Heater  voltage 

AV  =  ±  0.5  V 

Heater  current 

AI  =  ±  0.1  A 

Radial  acceleration 

Attr  =  ±  0.1-g 

Electrical  wire  resistance 

AiZwire  =  ±  2.4  X  10-3  n 

Calorimeter  inlet  temp. 

ATin  =  ±  0.07  K 

Calorimeter  outlet  temp. 

ATout  =  ±  0.08  K 

Evaporator  end  cap  temp. 

ATeec  =  ±  0.09  K 

Condenser  end  cap  temp. 

ATcec  =  ±  0.11  K 

Calculated  Values 

Heat  input 

See  Fig.  10 

Heat  transported 

AQt  =  ±  3.21  W 

Thermal  resistance 

See  Fig.  11 

Heat  transfer  coefficient 

See  Figs.  16  and  17 

All  of  the  measured  and  calculated  uncertainty  values  for  the  experimental  data  are 
presented  in  Table  12. 
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5  Anal3^ical  Capillary  Limit  Model  for  HGHPs 


The  analytical  model  used  in  the  present  work  was  proposed  by  Thomas  et  al. 
(1998).  The  model  used  an  iterative  approach  to  solve  for  the  heat  transported  by 
each  groove,  which  was  then  summed  to  determine  the  total  heat  transported  by  the 
heat  pipe.  Starting  with  the  following  equation  for  the  pressure  balance  in  a  heat 
pipe  (Faghri,  1995;  Chi,  1976) 


APcap,max  >  +  AP/  +  APbf 


an  equation  was  derived  for  the  pressure  balance  for  a  single  helical  groove  (Thomas 
et  al.,  1998) 


C  ^  (Le  +  2Z/ffl  +  Lc) 
J”c  hfg 


AfJiyQt  ^  t^lQg 
T^Pvrt  4pirfAgr(p 


-Pi 


{~5}6zi 


(28) 


An  extensive  analysis  of  the  geometric  parameters  of  the  helically-grooved  heat 
pipe  showed  that  the  groove  geometry  more  closely  matched  a  trapezoidal  groove  as 
opposed  to  the  rectangular  groove  assumed  by  Thomas  et  al.  (1998).  In  the  model 
by  Thomas  et  al.  (1998),  the  liquid  hydraulic  radius  corresponded  to  a  rectangular 
channel 


2hw 
2h  +  w 


(29) 


in  which  the  porosity  was  ip  =  1.0.  In  addition,  the  coefficient  of  friction  for  a 
rectangular  duct  proposed  by  Shah  and  Bhatti  (1987)  was  used 


fiRei  =  24(1  -  1.35530!  +  IMSla^  -  1.7012a®  +  0.9564a'‘  -  0.2537a®)  (30) 

where  the  aspect  ratio  is  a  =  w/h.  In  the  updated  model,  the  liquid  hydraulic  radius 
for  a  trapezoidal  groove  is 

_  2[wh  +  0.5h®(tan  di  +  tan  ^2)] 

~  w  +  h{sec  0i  +  sec  62) 
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The  porosity  is  defined  by  the  ratio  of  pore  volume  to  the  total  volume  of  the  wick 
structure.  Therefore,  for  a  trapezoidal  groove  the  porosity  is 

_Ys.  —  ^grLjAgr  _  Ngr[wh  +  0.5h^(tan 0i  +  tan  ^2)] 

^  Vt  Ltn{r‘^-r^)  7r(r2  -  r^) 

where  h,w,6i,  and  $2  are  shown  in  Fig.  1.  The  coefficient  of  friction  was  changed  to 

the  value  for  a  trapezoidal  duct  (/jRe;  =  14.21)  given  by  Shah  (1975). 

The  thermal  physical  properties  of  ethanol  given  by  Faghri  (1995)  were  used  by 
Thomas  et  al.  (1998).  It  was  found  from  an  extensive  literature  search  that  the 
properties  reported  by  Faghri  (1995)  were  for  aqueous  ethanol  which  is  96%  ethanol 
by  volume  and  4%  water.  In  the  present  experiment  99.9%  anhydrous  ethanol  was 
used  as  the  working  fluid.  Therefore,  the  density  of  the  liquid  and  vapor  will  be 
changed  to  the  best-fit  equations  found  in  Sec.  2.3. 

In  the  Thomas  et  al.  (1998)  model,  the  maximum  capillary  pressure  for  a  rectan¬ 
gular  groove  was  given  as 


cap, max 


where  a  is  the  surface  tension  and  Tc  is  the  capillary  pore  radius.  Thomas  et  al. 
(1999)  provides  the  capillary  pore  radius  for  liquid  flowing  in  trapezoidal  grooves.  As 
the  meniscus  recedes  into  the  groove  during  evaporation,  the  meniscus  bifurcates  into 
two  separate  flows  in  the  corners  of  the  groove.  This  point  of  bifurcation  is  assumed 
in  the  present  analysis  to  be  when  dryout  occurs  in  the  heat  pipe.  The  capillary  pore 
radius  at  the  point  of  bifurcation  is  (Thomas  et  al.,  1999) 

^  [1  +  (/5  +  tan  O)"^]  (34) 

where  P  =  w/2h  is  the  groove  aspect  ratio.  The  meniscus  contact  angle  at  bifurcation 
is  a  function  of  the  groove  geometry. 

tan-4  {^P^-^)cos9  +  p^cosSe  1 
[8;5cos0-l- (/52-|-4)sin^-t-/32sin30j 

Equations  (34)  and  (35)  were  used  to  determine  the  capillary  pore  radius  at  the  point 
of  dryout  in  the  following  manner.  The  groove  aspect  ratio  was  calculated  using  the 
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Table  13:  Comparison  of  the  parameters  used  in  the  analytical  model  for  predicting 
the  capillary  limit  of  the  helically-grooved  heat  pipe. 


Thomas  et  al.  (1998) 

Present 

h  (m) 

4.420x10-^ 

3.831  xl0"'‘ 

w  (m) 

4.735x10“^ 

3.445x10"'* 

La  (m) 

0.1016 

0.0823 

Tv  (m) 

6.730x10"'' 

6.793x10"^ 

P  (m) 

1.385 

1.358 

14.22 

14.21 

1.0 

0.5063 

n  (m) 

3.064x10"^ 

2.997x10"'* 

Th  (m) 

6.950x10"" 

6.985x10"*' 

Pv  (at  Tsat  =  40°C)  (kg/m^) 

0.2753 

0.3036 

Pi  (at  Tsat  =  40°C)  (kg/m^*) 

780.1 

772.3 

To  (m) 

4.735x10"'* 

2.252x10"'* 

groove  height  and  width  found  in  Table  3  (/?  =  0.4496).  As  per  the  discussion  by 
Thomas  et  al.  (1999),  this  value  of  ^  indicates  that  the  meniscus  must  recede  into  the 
groove  before  the  bifurcation  point  can  be  reached.  Therefore,  eqn.  (35)  was  solved 
iteratively  to  determine  the  aspect  ratio  for  which  the  meniscus  contact  angle  was 
equal  to  the  minimum  meniscus  contact  angle  for  the  copper-ethanol  system.  Faghri 
(1995)  gives  this  value  as  70  =  7°  for  a  receding  meniscus.  With  7^  =  70  =  T  and  6 
=  14.62°  (average  groove  half-angle  from  Table  3),  the  resulting  aspect  ratio  was  13 
=  1.2112.  Solving  for  h,  which  is  now  the  height  of  the  meniscus  above  the  groove 
bottom,  results  in  h  =  0.01422  cm.  Substituting  this  value  into  eqn.  (34)  for  the 
meniscus  radius  of  curvature  gives  Tg  =  0.0002252  cm. 

Table  13  shows  all  of  the  variables  used  by  Thomas  et  al.  (1998)  and  the  updated 
values  used  in  the  present  work.  The  most  significant  changes  are  to  the  wick  porosity 
and  capillary  pore  radius. 
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6  Results  and  Discussion 


The  objective  of  this  experiment  was  to  determine  the  steady-state  performance 
of  a  revolving  helically-grooved  heat  pipe  as  a  function  of  the  working  fluid  inventory. 
The  heat  input,  radial  acceleration  and  working  fluid  fill  were  varied  as  follows:  Qin 
=  10  to  180  W,  [url  =  0.01  to  10-g,  and  G  =  0.5,  1.0  and  1.5.  Thermocouples  on 
the  inboard,  outboard,  top,  and  bottom  sides  of  the  heat  pipe  (Fig.  9)  were  used  to 
determine  the  axial  and  circumferential  temperature  distributions.  Typical  steady- 
state  temperature  distributions  for  the  heat  pipe  for  G  =  1.0  at  |arl  =  0.01-g  are 
shown  in  Fig.  10.  All  of  the  steady-state  temperature  distributions  collected  during 
this  research  are  presented  in  Sec.  9.4.  For  low  power  input  levels,  the  temperature 
distribution  was  uniform.  As  the  power  input  increased,  the  temperatures  within 
the  evaporator  and  the  short  unheated  section  adjacent  to  the  evaporator  increased 
significantly,  indicating  a  partial  dryout  situation.  Since  the  coolant  temperature 
and  flow  rate  were  constant  for  all  tests,  the  adiabatic  and  condenser  temperatures 
increased  slightly  with  input  power.  Figure  11  shows  the  thermal  resistance  versus 
transported  heat  over  the  entire  range  of  radial  acceleration  for  each  fill  level.  In 
Fig.  11(a)  the  thermal  resistance  was  quite  high,  which  indicates  that  the  heat  pipe 
was  partially  dried  out  for  G  =  0.5,  even  at  the  lowest  power  input  levels.  However, 
the  thermal  resistance  decreased  significantly  as  the  radial  acceleration  increased, 
showing  that  the  capillary  pumping  ability  of  the  helical  grooves  increased.  For  G 
=  1.0  and  1.5,  the  thermal  resistance  decreased  and  then  increased  with  transported 
heat  when  dryout  commenced.  The  G  =  1.5  fill  tests  showed  dryout  occurring  only 
for  Icrl  =  0.01  and  2.0-g.  Dryout  was  not  reached  for  G  =  1.5  with  |or|  =  4.0,  6.0,  8.0 
and  10.0-g  due  to  reaching  the  maximum  allowable  heater  temperature  (Theater, max  = 
100°C).  The  capillary  limit  was  considered  to  be  reached  when  the  thermal  resistance 
began  to  increase  as  defined  by  Thomas  et  al.  (1998). 

Thomas  et  al.  (1998)  presented  a  mathematical  model  which  predicted  the  capil¬ 
lary  limit  of  a  helically-grooved  heat  pipe  subjected  to  a  transverse  body  force.  This 
model  accounted  for  the  geometry  of  the  heat  pipe  and  the  grooves  (including  helix 
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Rth  (K/W) 


Figure  12:  Capillary  limit,  versus  radial  acceleration  comparison  of  present  model  and 
Thomas  et  al.  (1998). 

pitch),  body  force  field  strength,  and  temperature-dependent  working  fiuid  proper¬ 
ties.  The  capillary  limit  versus  radial  acceleration  is  given  in  Fig.  12  for  various 
working  temperatures  with  the  Thomas  et  al.  (1998)  model  and  the  present  model  in 
which  improvements  were  made  as  discussed  in  Sec.  5.  The  capillary  limit  increased 
steadily  with  radial  acceleration  and  working  temperature.  The  present  model  shows 
a  significantly  lower  prediction  for  the  capillary  limit  when  compared  to  the  Thomas 
model  due  to  the  improved  geometric  measurements  and  working  fluid  property  equa¬ 
tions. 

Figure  13  shows  a  comparison  of  the  experimental  data  and  present  analytical 
model  for  the  capillary  limit  of  a  revolving  helically-grooved  heat  pipe.  No  attempt 
was  made  to  maintain  a  constant  adiabatic  temperature  during  the  experiments. 
Therefore,  the  working  fluid  temperature  in  the  model  was  set  to  the  adiabatic  tem¬ 
perature  found  experimentally.  For  G  =  0.5,  the  heat  pipe  operated  successfully  only 
for  \ar\  >  8.0-g.  In  Fig.  13(6)  {G  =  1.0),  the  capillary  limit  increased  significantly 
with  radial  acceleration.  With  the  heat  pipe  overfilled  by  50%  {G  =  1.5),  the  capillary 
limit  increased  dramatically,  showing  the  effect  that  overfilling  has  on  performance. 
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(m)  (M)  (m) 


Figure  13:  Comparison  of  present  model  and  experimental  capillary  limit  data  versus 
radial  acceleration:  (a)  G=  0.5;  (6)  G  =  1.0;  (c)  G  =  1.5. 
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The  agreement  of  the  analytical  model  was  very  good  for  G  =  1.0  as  expected.  For 
G  =  0.5,  the  model  overpredicted  the  experimental  data  because  it  was  assumed  that 
the  grooves  were  completely  filled.  For  G  =  1.5,  the  model  underpredicted  the  data 
due  to  the  assumption  that  no  liquid  communication  occurred  between  the  grooves. 

Temperatures  within  the  evaporator  section  are  shown  in  Figs.  14  and  15  for  |^|  = 
0.01-g  and  10.0-g,  respectively.  In  general,  the  temperatures  within  the  evaporator 
increased  with  transported  heat.  In  addition,  the  wall  temperatures  decreased  with 
G  for  a  given  heat  transport  due  to  the  fact  that  more  grooves  were  active.  The 
temperatures  along  the  length  of  the  evaporator  section  can  be  tracked  by  examining 
the  case  for  G  =  1.0.  Near  the  evaporator  end  cap,  the  temperatures  departed  those 
for  G  =  1.5  at  approximately  Qt  =  15  W  (Fig.  14(a)).  At  a;  =  92.1  mm  (Fig.  14(6)), 
this  departure  was  delayed  until  approximately  Qt  =  25  W,  and  at  a;  =  168  mm 
(Fig.  14(d)),  the  data  for  G  =  1.0  and  1.5  were  nearly  coincident.  This  behavior 
shows  that  the  grooves  were  essentially  full  near  the  adiabatic  section,  and  proceeded 
to  dry  out  closer  to  the  evaporator  end  cap,  as  expected.  Dryout  for  the  G  =  1.5 
case  can  be  seen  in  Fig.  14  (o)  where  the  temperatures  converged  to  nearly  the  same 
value  around  the  circumference.  It  should  be  noted  that  the  temperatures  around  the 
circumference  were  relatively  uniform  for  |ar|  =  0.01-g.  Evaporator  temperatures 
for  |dr|  =  10.0-g  are  shown  in  Fig.  15.  In  comparison  to  |dr|  =  0.01-g,  the  evaporator 
temperatures  were  in  general  lower  due  to  the  improved  pumping  ability  of  the  helical 
grooves  under  increased  radial  acceleration.  In  addition,  the  temperatures  tended  to 
overlap  over  a  greater  range  of  heat  transport  values.  In  contrast  to  |dr|  =  0.01-g, 
the  evaporator  temperature  variation  was  greater  around  the  circumference  at  higher 
Qt,  but  no  pattern  was  distinguishable  in  the  data. 

Local  heat  transfer  coefficient  data  versus  heat  transport  is  shown  in  Figs.  16 
and  17  for  \ar\  =  0.01-g  and  10.0-g.  Overall,  the  values  for  he  were  very  low  for  G 
=  0.5  due  to  the  fact  that  most  of  the  grooves  were  dried  out.  As  the  percent  fill 
increased  from  G  =  0.5  to  G  =  1.0,  the  heat  transfer  coefficient  increased  significantly. 
For  jorl  =  0.01-g  (Fig.  16),  hg  increased  and  then  decreased  with  transported  heat. 
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(Oo)  X  (Oo)  X  (Oo)  X  (Oo)  X 


(Do)  J,  (Oo)  X  (Oo)  X  (Oo)  X 


This  trend  was  also  reported  by  Vasiliev  et  al.  (1981)  for  an  aluminum  axially- 
grooved  heat  pipe  with  acetone  as  the  working  fluid.  For  (?  =  1.0  and  1.5,  the  heat 
transfer  coefficient  near  the  evaporator  end  cap  (Fig.  16(a))  decreased  until  all  of  the 
values  around  the  circumference  converged.  Closer  to  the  adiabatic  section,  the  heat 
transfer  coefficient  values  around  the  circumference  had  not  yet  converged,  showing 
these  portions  to  still  be  active.  For  |5r|  =  10.0-g  (Fig.  17),  the  values  of  were 
significantly  more  uniform  around  the  circumference  and  along  the  axial  direction, 
even  during  a  dryout  event  (G  =  1.0,  Fig.  17(o)).  In  addition,  the  heat  transfer 
coefficient  seems  to  be  more  constant  with  respect  to  the  transported  heat  compared 
to  jcrl  =  0.01-g.  During  the  experiments,  the  heat  pipe  working  temperature  was 
not  constant,  which  resulted  in  changes  in  the  specific  volume  of  the  liquid  and  vapor 
of  the  working  fluid.  Since  the  heat  pipe  was  filled  at  room  temperature,  it  was 
important  to  quantify  the  potential  effects  of  the  change  in  volume  of  liquid  in  the 
grooves  with  temperature.  Figure  18  shows  the  variation  of  the  percentage  of  groove 
volume  occupied  by  liquid  G  with  saturation  temperature  for  the  three  fill  amounts 
over  the  range  of  working  temperatures  seen  in  the  experiments.  The  maximum 
percent  difference  was  2.7%,  which  was  not  deemed  to  be  significant. 
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(W/m^K)  he  (W/m^K)  he  (W/m^K)  he  (W/m^K) 


168  mm. 


7  Conclusions 


The  effect  of  fluid  inventory  on  the  steady-state  performance  of  a  helically-grooved 
copper-ethanol  heat  pipe  has  been  examined  both  experimentally  and  analytically. 
It  was  found  that  the  capillary  limit  increased  and  the  thermal  resistance  decreased 
significantly  as  the  amount  of  working  fluid  within  the  heat  pipe  increased.  In  ad¬ 
dition,  the  evaporative  heat  transfer  coeflScient  was  found  to  be  a  strong  function 
of  the  fill  amount.  The  updated  analytical  model  was  in  very  good  agreement  with 
the  experimental  capillary  limit  results  for  (?  =  1.0.  However,  the  analytical  model 
overpredicted  the  capillary  limit  data  for  G  =  0.5  and  underpredicted  the  data  for  G 
=  1.5. 
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9  Appendices 


9.1  Heat  Pipe  Filling  Station  Operating  Procedure 

In  this  section,  the  procedure  for  filling  a  heat  pipe  is  outlined.  Table  14  shows 
the  valving  sequence  during  a  typical  filling  procedure.  Prior  to  filling,  the  heat  pipe 
container  was  evacuated  for  24  hr  using  a  turbomolecular  vacuum  pump  to  a  pressure 
of  10“®  Torr.  The  evacuated  heat  pipe  was  then  installed  onto  the  filling  station  at 
valve  5. 

1.  Close  all  valves. 

2.  Fill  the  working  fluid  reservoir  with  500  ml  of  fluid  and  remove  any  noncon¬ 
densable  gas  from  the  fluid  by  the  freeze/thaw  method.  Prior  to  handling  LN2, 
read  the  Material  Safety  Data  Sheet  (MSDS),  and  wear  protective  clothing  and 
goggles. 

(a)  Fill  the  LN2  dewax.  Connect  a  vacuum  line  from  the  manifold  drain  plug 
port  to  valve  3. 

(b)  Slowly  raise  the  LN2  dewar  with  the  lab  jack  until  the  working  fluid  reser¬ 
voir  is  partially  submerged  to  freeze  the  working  fluid.  The  pressure  read 
by  the  vacuum  gage  will  be  less  than  100  jumHg  when  the  ethanol  freezes. 

(c)  Open  and  then  close  valves  3  and  7  for  10  min.  to  clear  the  reservoir  of 
noncondensable  gases.  Lower  the  LN2  dewar  to  thaw  the  working  fluid. 

(d)  Repeat  steps  2(6)  and  2(c)  at  least  three  times. 

3.  Remove  the  air  from  the  fill  station  manifold  and  heat  pipe  fill  valve  stem. 

(a)  Install  the  manifold  drain  plug  and  the  evacuated  heat  pipe  on  the  fill 
station. 

(b)  Open  valves  2,  5,  6  and  7.  Evacuate  the  system  for  20  minutes  (Manifold 
pressure  <  60  /im). 
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4.  Fill  the  manifold  and  dispensing  burette  completely  with  fluid. 

(a)  Isolate  the  manifold  from  the  vacuum  source  by  closing  valves  6  and  7. 

(b)  Slowly  open  valve  4  to  draw  fluid  into  the  dispensing  burette. 

(c)  If  necessary,  pressurize  the  fluid  reservoir  by  attaching  the  GN2  feed  line 
to  valve  3. 

(d)  Open  the  valve  on  the  GN2  cylinder  and  adjust  the  regulated  pressure  of 
the  feed  line  to  less  than  10  psig. 

(e)  Open  valve  3  and  slowly  open  and  close  valve  1  to  completely  fill  the 
dispensing  burette  with  fluid.  Vapor  above  the  fluid  will  be  vented  to 
atmosphere. 

(f)  First  close  valve  4,  then  close  valve  3  and  the  valve  on  the  GN2  cylinder. 

5.  Purge  the  manifold  of  vapor  bubbles  by  pressurizing  the  fluid  in  the  dispensing 
burette  using  the  GN2  feed  line  attached  to  the  pressure  gage  (dashed  line  in 
Fig.  6). 

(a)  Open  the  valve  on  the  GN2  cylinder  and  adjust  the  pressure  to  less  than 
5  psig  (GN2  will  escape  from  the  1  psig  pressure  relief  valve).  Open  valve 
1  and  remove  the  manifold  drain  plug. 

(b)  Open  and  then  close  valve  6  to  drain  down  the  fluid  to  the  first  graduation 
mark  on  the  dispensing  burette.  Fluid  will  appear  at  the  manifold  drain. 

(c)  Close  valve  2  and  note  the  change  in  the  fluid  meniscus  level  when  opening 
and  closing  the  pressure  relief  valve. 

(d)  Open  valve  2  and  again  note  the  change  in  the  meniscus  level  when  opening 
and  closing  the  pressure  relief  valve.  If  there  is  a  significant  difference 
between  the  change  found  in  steps  5(c)  and  5(d),  the  presence  of  trapped 
vapor  bubbles  in  the  manifold  is  indicated. 
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(e)  Close  and  open  valves  2  and  5  several  times  while  cycling  the  the  pressure 
relief  valve  to  remove  any  trapped  vapor  bubbles  within  these  valves.  The 
volume  displaced  by  each  valve  (without  trapped  vapor)  is  0.06  ml.  Valves 
2  and  5  must  be  open  at  the  end  of  this  sub-process. 

6.  Dispense  the  working  fluid  into  the  heat  pipe. 

(a)  Record  the  dispensing  burette  meniscus  level. 

(b)  Slowly  open  and  then  close  valve  8  (heat  pipe  fill  valve)  to  meter  in  the  de¬ 
sired  volume  of  fluid  into  the  heat  pipe  by  noting  the  drop  in  the  meniscus 
level. 

(c)  Again  record  the  dispensing  burette  meniscus  level. 

(d)  Close  valve  5  and  remove  the  filled  heat  pipe. 

7.  Purge  the  filling  station  of  working  fluid. 

(a)  Open  valve  3  to  atmosphere  and  then  slowly  open  valve  4  to  purge  the 
unused  fluid  in  the  dispensing  burette  back  into  the  fluid  reservoir. 

(b)  Close  valve  4  and  then  valve  3,  and  remove  the  working  fluid  reservoir  from 
the  filling  station  at  valve  4. 

(c)  Open  valves  4,  5  and  6  to  allow  the  GN2  to  evaporate  any  remaining 
working  fluid  within  the  filling  station. 

(d)  Once  the  manifold  is  dry,  close  the  valve  on  the  GN2  cylinder,  then  close 
valves  1,  4,  5,  and  6  and  install  the  manifold  drain  plug. 
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Table  14:  Flow  chart  for  fill  station  operating  procedure  (0  =  Open,  C  =  Closed, 
OC  =  Open  then  Close,  CO  =  Close  then  Open). 


Step 

1 

2 

3  1 

4  1  5  1 

Valve 

6  1  7  1  8  1  M.D.P.  1  GN2  I  P.R.V. 
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C 

C 

m 

c 

c 

0 

c 
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9.2  Thermocouple  Calibration  Procedure 

The  thermocouples  used  in  the  present  experiment  were  calibrated  from  20  to 
100®C  in  increments  of  10°C  with  a  Hart  Scientific  (Model  5691)  NIST-Traceable 
resistance  temperature  detector  (RTD)  and  a  low-dead  band  (±  0.01°C)  Brinkmann 
(Model  RC20)  recirculating  bath.  The  thermocouples  were  mounted  to  the  RTD 
probe  and  submerged  into  the  recirculating  bath.  The  recirculation  of  the  bath 
oil  (Brayco  micronic  889)  caused  significant  temperature  gradients  throughout  the 
bath.  To  address  this  problem  a  small  can  was  submerged  in  the  bath  and  the 
RTD/thermocouple  bundle  was  placed  inside  the  can.  Electrical  noise  was  found 
to  be  present  in  the  test  cell  where  the  thermocouples  were  calibrated,  as  shown 
in  Fig.  19(a).  To  address  this  problem  the  data  was  averaged  every  10  s  over  the 
entire  sampling  interval.  Figure  19(6)  shows  a  typical  averaged  data  set  taken  at  a 
setpoint  temperature  of  90°C  for  3.5  hr,  where  the  uncertainty  of  the  readings  has 
been  significantly  reduced.  Even  though  the  averaged  data  was  relatively  constant, 
the  recirculating  bath  was  stabilized  for  2  hr  before  data  was  taken  at  each  setpoint 
temperature  to  ensure  that  steady-state  conditions  had  been  reached. 

Since  some  of  the  thermocouples  were  switched  by  relays,  a  test  was  performed  in 
order  to  ensure  that  switching  the  relays  on  and  off  did  not  affect  the  thermocouple 
readings.  The  test  consisted  of  starting  the  data  acquisition  and  waiting  1  min., 
switching  the  relays  to  the  on  position  and  waiting  15  min.,  then  switching  the  relays 
to  the  off  position  and  waiting  15  min.  The  result  of  the  test  showed  there  was  no 
significant  change  in  the  temperature  with  time  for  any  of  the  thermocouples.  A 
series  of  tests  was  also  performed  on  the  thermocouple  relays  to  determine  the  effect 
the  relay  contacts  had  on  the  thermocouple  output  temperature.  Some  of  the  tests 
consisted  of  blowing  hot  and  cold  air  across  the  relay  contacts.  The  results  of  these 
tests  showed  that  the  relays  had  to  be  insulated  and  covered  with  a  box  to  keep  the 
convective  heat  losses  to  a  minimum  while  the  centrifuge  table  was  spinning. 

The  following  is  the  procedure  used  in  the  calibration  process. 
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Temperature  by  thermocouple 


1.  Prepare  the  thermocouple  coils  to  match  the  conditions  during  testing  by  en¬ 
suring  they  are  insulated  and  all  covers  are  in  place. 

2.  Mount  all  of  the  thermocouples  (including  the  calorimeter  thermocouples)  onto 
the  RTD  probe  as  close  to  the  sensing  portion  as  possible. 

3.  Place  the  RTD/thermocouple  bundle  vertically  into  the  recirculating  bath  with 
a  laboratory  stand.  Submerge  the  bundle  as  far  as  possible  without  the  RTD 
or  any  of  the  thermocouples  touching  the  bottom  or  sides  of  the  recirculating 
bath. 

4.  Turn  on  the  recirculating  bath  and  adjust  the  setpoint  to  20“C. 

5.  Turn  on  the  centrifuge  table  power. 

6.  Turn  on  the  RTD  computer  and  initiate  the  program  associated  with  the  RTD. 

7.  Turn  on  the  control  room  electronics  and  thermocouple  data  acquisition  pro¬ 
gram  Viewdao. 

8.  Place  the  thermocouple  coil  switch  in  the  OFF  position. 

9.  Open  Viewdac  sequences  calibrat.beq  and  calibrwl.beq. 

10.  Use  calibrat.beq  to  monitor  the  bath  temperature  and  calibrwl.beq  to  write 
data  to  a  file. 

11.  Ensure  that  the  sample  rate  is  set  to  10  Hz  and  the  output  file  name  has  proper 
identification  for  the  thermocouple  coil  switch  in  the  OFF  position. 

12.  Start  the  data  acquisition  loop  and  run  for  120  s. 

13.  Record  an  RTD  value  manually  every  ten  seconds  while  the  data  acquisition 
loop  is  running. 

14.  Change  the  output  file  name  to  denote  that  the  thermocouple  coil  switch  is  in 
the  ON  position. 
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15.  Switch  the  thermocouple  coil  switch  to  the  ON  position.  Wait  ten  seconds  to 
ensure  coil  stability. 

16.  Start  the  data  acquisition  loop  and  run  for  120  s. 

17.  Record  an  RTD  value  manually  every  ten  seconds  while  the  data  acquisition 
loop  is  running. 

18.  Switch  the  thermocouple  coil  switch  to  the  OFF  position. 

19.  Advance  the  recirculating  bath  to  the  next  setpoint  temperature. 

20.  Repeat  steps  12-19  for  the  next  setpoint. 

Coefficients  for  the  linear  curve  fit  for  each  thermocouple  are  given  in  Table  15 
for  the  range  20°  C  ^  Treading  ^  100°. 

TjwAual  ~  “b  O^lTejiding  ^  ^ 

Using  the  standard  error  given  by  Coleman  and  Steele  (1989),  the  uncertainty  asso¬ 
ciated  with  each  of  the  thermocouple  calibration  equations  is 
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Table  15:  Coefficients  for  thermocouple  calibration. 


Thermocouple 

ao 

Oil 

€ 

Data  Acquisition 
Channel 

Coolant  Inlet  (11) 

0.0699 

TC3 

Coolant  Outlet  (12) 

-0.8748547 

0.9996509 

0.0841 

TC4 

1-Inboard 

-0.0521081 

0.9959866 

0.0732 

TCO 

2-Ihboard 

-0.2641292 

1.0001844 

0.0800 

TC5 

3-Inboard 

-0.8430881 

1.0003875 

0.0760 

TC6 

4-Inboard 

-0.5082865 

0.9970537 

0.0736 

TC7 

5-Inboard 

-0.3317360 

0.9991215 

0.0702 

TC8 

6-Inboard 

-0.5198530 

0.9961772 

0.0828 

TC13 

8-Inboard 

-0.3119424 

0.9961133 

0.0732 

TC14 

10-Inboard 

-0.6501503 

0.9965431 

0.0686 

TC2 

,  1-Outboard 

-0.5721585 

0.9959817 

0.0776 

TCI 

2-Outboard 

-0.1469691 

0.9960199 

0.0800 

TC5 

3-Outboard 

-0.8401706 

1.0008819 

0.0769 

TC6 

4-6utboard 

-0.4942304 

0.9970113 

0.0699 

TC7 

5-Outboard 

-0.3201927 

0.9992546 

0.0748 

TC8 

7-Outboard 

-0.5018720 

0.9960913 

0.0816 

TC13 

9-Outboard 

-0.2943188 

0.9959160 

0.0750 1 

TC14 

lO-Outboard 

0.0738 

TC15 

2-Top 

-0.5137903 

1.0005327 

0.1341 

TCO 

3-Top 

-0.2518604 

1.0001036 

0.0759 

TCIO 

4-Top 

-0.3904948 

1.0013209 

0.0656 

TCll 

5-Top 

0.4147904 

0.9950726 

0.0649 

TC12 

2-Bottom 

-0.5077360 

1.0004865 

0.0661 

TCO 

3-Bottom 

-0.2490069 

TCIO 

4-Bottom 

-0.4025206 

1.0014220 

0.0686 

TCll 

5-Bottom 

0.0652 

TC12 
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9.3  Calculated  Experimented  data 
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Table  16:  Evaporative  heat  transfer  coefficient  (W/m^K)  for  thermocouple  2,  G  = 
0.5  and  =  0.01-g. 


Bottom 

Outboard 

Top 

Inboard 

Qin  (W) 

Hq 

Ahe 

/ig 

/ig 

<1 

/ig 

Ahe 

11.254 

54.580 

53.971 

28.184 

54.737 

28.577 

15.494 

39.413 

15.219 

39.808 

15.374 

38.956 

15.044 

39.832 

15.382 

20.777 

26.511 

9.834 

26.141 

9.698 

26.821 

9.950 

25.045 

28.554 

7.966 

28.740 

8.018 

28.158 

7.856 

28.896 

8.062 

Table  17:  Evaporative  heat  transfer  coefficient  (W/m^K)  for  thermocouple  B,  G  = 
0.5  and  ar  =  0.01-g. 


Bottom 

Outboard 

Top 

Inboard 

Qin  (W) 

/ig 

i^n 

/ig 

Ahe 

/ig 

■aBI 

11.254 

71.586 

37.373 

34.826 

66.082 

34.500 

15.494 

44.767 

16.502 

42.918 

16.575 

20.777 

28.368 

27.366 

10.152 

27.679 

10.268 

25.045 

29.829 

8.322 

29.567 

8.249 

28.894 

8.061 

29.349 

8.188 

Table  18:  Evaporative  heat  transfer  coefficient  (W/m^K)  for  thermocouple  4,  G  = 
0.5  and  ar  =  0.01-g. 


Bottom 

Outboard 

Top 

Inboard 

Qln  (W) 

mam 

Ahe 

/ig 

mm 

/ig 

mm 

/ig 

11.254 

238.999 

125.195 

184.815 

96.720 

165.971 

86.843 

178.979 

93.557 

24.907 

66.693 

25.761 

64.852 

25.047 

66.219 

25.576 

20.777 

34.800 

12.910 

35.683 

12.996 

35.804 

13.282 

25.045 

34.787 

9.705 

35.581 

9.926 

35.051 

9.778 

35.799 

9.987 
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Table  19:  Evaporative  heat  transfer  coefficient  (W/m^K)  for  thermocouple  5,  G  = 
0.5  and  Or  =  0.01-g. 


Bottom 

Outboard 

Top 

Inboard 

Qin 

mm 

mm 

/ig 

mam 

/ig 

<1 

660.165 

611.749 

328.146 

788.292 

423.748 

920.380 

501.517 

15.494 

876.782 

350.646 

400.343 

155.831 

359.535 

139.983 

53.186 

92.298 

34.257 

90.234 

33.488 

25.045 

90.659 

25.309 

72.417 

20.210 

72.585 

20.254 

69.268 

19.328 

Table  20:  Evaporative  heat  transfer  coefficient  (W/m^K)  for  thermocouple  2,  G  = 
1.0  and  ttr  =  0.01-g. 


Bottom 

Outboard 

Top 

Inboard 

Ahe 

/ig 

/irg 

Ahe 

/Zrg 

Ahe 

10.903 

557.26 

257.17 

568.41 

265.15 

21.176 

622.95 

124.73 

608.62 

122.12 

498.53 

100.61 

532.93 

107.02 

30.494 

277.81 

37.22 

290.38 

283.08 

38.04 

40.053 

142.35 

14.75 

146.50 

15.19 

145.46 

15.20 

49.903 

100.06 

9.00 

98.77 

8.89 

100.47 

9.04 

Table  21:  Evaporative  heat  transfer  coefficient  (W/m^K)  for  thermocouple  3,  G  = 
1.0  and  Or  =  0.01-g. 


Bottom 

Outboard 

Top 

Inboard 

Qin  (W) 

Hq 

Ahe 

Hq 

Ahe 

/ig 

Ahe 

/ig 

Ahe 

10.903 

1402.98 

718.49 

704.90 

525.04 

955.67 

468.29 

21.176 

1674.12 

382.68 

816.53 

166.82 

1124.19 

236.99 

1127.08 

236.72 

30.494 

1325.08 

197.80 

676.35 

94.98 

mm 

40.053 

397.60 

42.54 

301.54 

31.38 

263.87 

27.41 

49.903 

149.08 

13.43 

147.24 

13.26 

138.36 

12.45 

138.49 

12.46 
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Table  22:  Evaporative  heat  transfer  coeflBcient  (W/m^K)  for  thermocouple  4,  G  = 
1.0  and  Or  =  0.01-g. 


Qin  (W) 


10.903 


21.176 


30.494 


40.053 


49.903 


Bottom 


1095.94 


1347.75 


1354.52 


1248.87 


847.62 


Ah, 


542.02 


287.24 


207.02 


146.83 


79.50 


Outboard 


ho 


911.73 


1131.17 


1038.21 


827.96 


472.47 


Ah, 


444.91 


241.54 


152.70 


90.11 


43.12 


Top 


h. 


1018.26 


1325.19 


1315.68 


987.30 


477.35 


Ah, 


496.46 


289.71 


192.11 


106.86 


43.42 


Inboard 


h. 


1499.84 


1634.73 


2015.83 


1588.34 


687.18 


Ah, 


826.94 


383.83 


337.37 


194.26 


64.25 


Table  23:  Evaporative  heat  transfer  coefficient  (W/m^K)  for  thermocouple  5,  G  = 
1.0  and  Or  =  0.01-g. 


Bottom 

Outboard 

Top 

Inboard 

Qiu 

Ah, 

/ig 

Ah, 

Hq 

Ah, 

/ig 

Ah, 

10.903 

1248.05 

618.59 

785.68 

374.66 

894.79 

423.76 

1174.16 

578.17 

21.176 

1239.50 

271.27 

928.74 

191.37 

1035.20 

220.32 

1191.51 

255.87 

30.494 

1265.95 

184.38 

921.92 

128.64 

1083.38 

154.07 

1204.23 

174.67 

40.053 

1192.77 

132.43 

887.30 

97.74 

1043.66 

114.21 

1140.45 

127.81 

49.903 

1122.63 

111.39 

849.51 

80.00 

980.95 

93.10 

1076.25 

103.55 

Table  24:  Evaporative  heat  transfer  coefficient  (W/m^K)  for  thermocouple  2,  G  = 
1.5  and  Ur  =  0.01-g. 


Bottom 

Outboard 

Top 

Inboard 

Qin  (W) 

/ig 

Z^/ig 

/ig 

<1 

/ig 

Ah, 

Ah, 

10.124 

987.10 

368.44 

775.69 

292.30 

898.09 

333.98 

1051.70 

199.72 

29.462 

1111.74 

93.65 

953.77 

117.78 

39.808 

797.77 

74.96 

801.47 

75.54 

687.69 

65.10 

50.142 

794.32 

63.26 

688.77 

54.91 

60.577 

652.89 

40.23 

629.16 

30.22 

579.27 

35.67 

584.59 

31.87 

557.571 

521.71 

28.12 

80.619 

532.34“1 

24.67 

502.15 

476.34 

21.84 

89.645 

540.57 

31.18 

511.88 

484.42 

23.32 

101.692 

438.04 

18.83 

405.92 

16.73 

110.802 

307.01 

11.75 

316.84 

12.21 

307.40 

12.03 

312.52 

12.19 

Table  25:  Evaporative  heat  transfer  coefficient  (W/m^K)  for  thermocouple  3,  G  = 
1.5  and  =  0.01-g. 


Bottom 

Outboard 

Top 

Inboard 

«»  (W) 

he 

Hq 

Ahe 

Hq 

/ig 

WEIM\ 

10.124 

1569.96 

650.12 

964.54 

365.12 

1262.58 

497.12 

836.24 

307.83 

20.024 

1730.98 

364.55 

978.85 

185.45 

1460.90 

188.77 

29.462 

1904.73 

271.33 

1039.95 

130.69 

1500.52 

145.53 

141.25 

843.85 

80.73 

1020.95 

98.50 

819.04 

77.93 

50.142 

1071.25 

89.56 

663.19 

52.30 

782.14 

61.31 

656.62 

52.28 

60.577 

714.24 

52.92 

526.12 

32.67 

592.79 

37.20 

513.79 

32.27 

69.754 

615.53 

38.50 

473.85 

25.41 

528.96 

28.88 

464.91 

26.22 

80.619 

549.94 

28.33 

435.76 

19.89 

485.00 

22.76 

429.34 

20.16 

89.645 

673.87 

45.70 

498.63 

25.10 

557.01 

27.56 

499.08 

24.48 

101.692 

645.42 

31.70 

488.83 

19.92 

539.80 

21.97 

481.49 

mm 

110.802 

671.19 

26.74 

488.41 

18.63 

526.78 

20.82 

479.25 

18.37 

Table  26:  Evaporative  heat  transfer  coefficient  (W/m^K)  for  thermocouple  A,  G  = 
1.5  and  =  0.01-g. 


Bottom 

Outboard 

Top 

Inboard 

«i.  (W) 

he 

Ahe 

Ahe 

Hq 

Ahe 

<1 

10.124 

1224.23 

462.42 

1264.21 

497.76 

1264.24 

666.10 

20.024 

1293.35 

251.80 

1230.76 

240.52 

1425.45 

288.04 

1705.03 

362.79 

29.462 

1399.15 

182.91 

1181.53 

154.20 

1489.83 

260.32 

39.808 

1361.70 

133.50 

1108.98 

111.61 

1416.39 

138.90 

1906.83 

213.24 

50.142 

1054.00 

86.27 

909.67 

74.72 

1037.44 

83.96 

1438.24 

127.36 

60.577 

789.57 

50.43 

728.30 

45.64 

720.66 

45.25 

855.04 

55.53 

69.754 

637.53 

34.65 

617.56 

33.72 

570.07 

30.73 

607.55 

33.45 

80.619 

538.84 

24.57 

532.37 

24.19 

480.59 

21.84 

504.22 

23.15 

89.645 

592.38 

34.37 

566.40 

31.23 

523.88 

24.94 

560.89 

26.98 

101.692 

644.04 

582.55 

23.70 

550.08 

22.21 

620.25 

25.05 

110.802 

655.10 

25.51 

575.90 

22.24 

552.01 

21.55 

646.52 

25.14 
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Table  27:  Evaporative  heat  transfer  coefficient  (W/m^K)  for  thermocouple  5,  G  = 
1.5  and  Or  =  0.01-g. 


Table  28:  Evaporative  heat  transfer  coefficient  (W/m^K)  for  thermocouple  2,  G  = 
0.5  and  ar  =  10.0-g. 
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Table  29:  Evaporative  heat  transfer  coefficient  (W/m^K)  for  thermocouple  3,  G  = 
0.5  and  =  10.0-g. 


Bottom 

Outboard 

Top 

Inboard 

Hq 

MSM 

mm 

10.735 

729.30 

177.96 

579.12 

139.84 

594.26 

143.66 

15.694 

70.54 

467.42 

71.60 

337.05 

51.15 

32.90 

270.51 

32.63 

256.09 

30.85 

211.73 

25.44 

14.17 

122.51 

12.74 

29.973 

95.08 

10.21 

93.22 

10.01 

89.24 

9.59 

89.08 

9.57 

33.634 

79.63 

7.36 

78.11 

7.22 

75.60 

6.99 

75.95 

7.02 

40.170 

64.03 

5.09 

63.25 

5.03 

61.24 

4.87 

61.80 

4.91 

43.285 

60.19 

4.67 

59.76 

4.64 

57.84 

4.49 

58.38 

4.53 

Table  30:  Evaporative  heat  transfer  coefficient  (W/m^K)  for  thermocouple  4,  G  — 
0.5  and  =  10.0-g. 


Bottom 

Outboard 

Top 

Inboard 

ein  (W) 

/ig 

Ahe 

/ig 

Ahe 

Hq 

Ahe 

10.735 

1009.38 

254.97 

774.51 

190.40 

663.68 

177.73 

15.694 

927.33 

148.84 

746.73 

117.68 

627.76 

96.92 

20.869 

597.21 

73.70 

559.62 

69.15 

472.33 

57.68 

25.298 

285.69 

34.61 

276.29 

27.18 

29.973 

217.31 

23.44 

189.07 

20.34 

33.634 

148.04 

13.72 

163.74 

15.20 

146.44 

13.56 

136.30 

12.61 

40.170 

104.50 

8.31 

110.33 

8.78 

104.12 

8.28 

99.23 

7.89 

43.285 

93.07 

7.23 

96.69 

7.52 

93.44 

7.26 

89.86 

6.98 
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Table  31:  Evaporative  heat  transfer  coefficient  (W/m^K)  for  thermocouple  5,  G  = 
0.5  and  Cr  =  10.0-g. 


Table  32:  Evaporative  heat  transfer  coefficient  (W/m^K)  for  thermocouple  2,  G  = 
1.0  and  Of  =  10.0-g. 
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Table  33:  Evaporative  heat  transfer  coeflBcient  (W/m^K)  for  thermocouple  3,  G  = 
1.0  and  Cr  =  10.0-g. 


Table  34:  Evaporative  heat  transfer  coefficient  (W/m^K)  for  thermocouple  4,  G  = 
1.0  and  Cr  =  10.0-g. 
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Table  35:  Evaporative  heat  transfer  coefficient  (W/m^K)  for  thermocouple  5,  G 
1.0  and  ttr  =  10.0-g. 


Bottom 


Outboard 


Top 


Inboard 


Q 


in 


he 


Ahe 


he 


Ahe 


Ahe 


he 


Ahe 


10.966 


717.21 


185.83 


903.80 


240.71 


758.06 


195.90 


960.79 


255.12 


20.617 


1114.11 


145.03 


648.33 


81.44 


721.99 


90.66 


909.16 


116.66 


29.937 


1148.33 


102.74 


652.49 


55.40 


708.48 


60.22 


869.98 


76.01 


40.167 


1159.28 


76.00 


652.43 


40.91 


698.94 


44.10 


843.92 


54.13 


49.244 


1162.59 


62.35 


649.35 


33.29 


683.72 


34.92 


836.54 


43.42 


60.373 


1158.07 


53.23 


648.36 


28.83 


674.15 


29.70 


811.23 


36.11 


70.243 


1152.33 


46.91 


654.46 


25.62 


676.02 


26.52 


806.42 


31.92 


80.156 


1129.35 


42.24 


653.01 


23.58 


673.65 


24.27 


798.64 


29.13 


88.516 


1155.44 


38.64 


649.23 


21.10 


665.69 


21.72 


787.13 


26.02 


98.323 


1127.79 


36.04 


639.64 


20.05 


658.06 


20.53 


771.64 


24.24 


111.223 


1132.65 


34.66 


649.71 


19.18 


669.96 


19.76 


780.17 


23.13 


119.747 


1119.88 


32.42 


646.48 


667.68 


18.82 


775.49 


131.622 


1132.93 


31.58 


651.00 


676.16 


18.29 


781.98 


139.843 


1087.71 


30.19 


627.91 


647.74 


17.46 


746.76 


22.09 


21.37 


20.26 


18.20 

17.62 

16.95 


Table  36:  Evaporative  heat  transfer  coefficient  (W/m^K)  for  thermocouple  2,  G  = 
1.5  and  =  10.0-g. 


Bottom 

Outboard 

Top 

Inboard 

<?i.  (W) 

Ahe 

Hq 

/ig 

<1 

/ig 

KB 

10.281 

449.88 

107.45 

507.44 

122.39 

395.13 

95.90 

418.71 

100.39 

20.038 

421.66 

51.60 

477.88 

58.98 

380.43 

50.35 

30.875 

397.20 

35.02 

365.54 

28.53 

395.26 

30.59 

41.022 

404.27 

55.40 

463.62 

63.79 

366.42 

68.49 

395.65 

73.51 

50.181 

427.01 

49.21 

360.19 

60.10 

61.273 

33.50 

361.90 

38.97 

25.71 

448.51 

29.59 

371.39 

32.97 

407.34 

35.80 

80.548 

394.49 

13.79 

454.15 

15.94 

368.87 

14.21 

89.981 

394.57 

12.92 

15.07 

369.10 

13.33 

100.278 

403.17 

12.16 

465.86 

14.16 

379.79 

12.64 

111.715 

435.63 

19.55 

512.23 

22.39 

417.47 

18.18 

455.13 

19.71 

122.090 

462.96 

19.67 

555.08 

23.01 

454.68 

18.63 

482.78 

19.77 

129.122 

472.64 

12.61 

571.10 

15.31 

464.73 

13.18 

mma^ 

15.50 

480.83 

13.48 

149.761 

12.64 

602.03 

15.28 

481.68 

12.29 

516.89 

13.13 

■lianrii 

16.24 

500.24 

iiggiia 

13.36 

16.29 

511.78 

12.88 

12.98 

660.12 

16.09 

511.08 

13.30 
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Table  37:  Evaporative  heat  transfer  coefficient  (W/m^K)  for  thermocouple  Z,  G  — 
1.5  and  Ur  =  10.0-g. 


Bottom 

Outboard 

Top 

Inbo 

ard 

Ob.(W) 

Ahe 

Hq 

Ahe 

/ig 

Ahe 

/ig 

Ahe 

10.281 

577.30 

139.40 

446.53 

107.27 

507.09 

121.81 

393.17 

94.08 

20.038 

491.65 

60.41 

414.94 

50.89 

475.89 

58.58 

376.39 

45.97 

30.875 

422.00 

28.60 

428.94 

33.24 

341.93 

26.41 

41.022 

445.48 

61.17 

392.22 

53.73 

434.92 

85.17 

350.74 

66.30 

50.181 

395.84 

45.55 

354.53 

40.71 

418.53 

67.31 

337.12 

52.21 

61.273 

390.18 

28.85 

354.22 

26.14 

396.00 

40.26 

322.65 

32.04 

70.586 

378.91 

24.86 

346.81 

36.97 

326.00 

29.08 

80.548 

377.81 

13.20 

350.07 

13.92 

320.21 

11.18 

89.981 

365.13 

12.02 

341.24 

11.20 

388.39 

12.68  ^ 

311.81 

10.17 

100.278 

426.12 

12.94 

406.60 

12.60 

456.22 

14.02 

354.11 

10.73 

111.715^ 

467.72 

23.57^ 

454.32 

24.90 

520.15 

22.58 

387.43 

16.78 

122.090 

478.20l 

23.04 

454.51“] 

23.67" 

541.05 

22.20 

393.12 

16.04 

129.122 

506.96 

14.00 

482.42 

13.34 

592.51 

16.56 

413.34 

11.22 

140.716 

536.94 

14.28 

619.94 

16.79 

434.71 

11.58 

149.761 

555.63 

14.38 

543.26 

13.77 

656.71 

16.82 

446.01 

11.26 

158.110 

588.75 

15.59 

555.49 

14.28 

679.36 

18.12 

461.34 

11.94 

169.488 

603.73 

15.11 

590.11 

15.17 

726.39 

18.85 

472.10 

11.76 

180.974 

620.00 

15.07 

619.34 

14.99 

750.54 

19.06 

479.35 

11.80 
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Table  38:  Evaporative  heat  transfer  coefficient  (W/m^K)  for  thernaocouple  4,  G 
1.5  and  =  10.0-g. 


Bottom 

Outboard 

Top 

Inboard 

<3to(W) 

/ig 

/ig 

Ahe 

/ig 

WEM 

he 

Ah, 

10.281 

928.13 

231.36 

639.44 

156.02 

621.05 

150.31 

797.22 

198.45 

20.038 

918.59 

117.63 

629.37 

78.24 

596.45 

30.875 

748.58 

59.57 

550.11 

43.12 

513.69 

39.96 

567.21 

44.60 

41.022 

759.48 

561.64 

77.73 

541.74 

122.58 

50.181 

606.42 

70.64 

490.18 

56.65 

470.67 

87.62 

61.273 

530.68 

39.82 

467.79 

34.69 

47.34 

32.45 

450.13 

29.67 

426.14 

40.53 

432.28 

42.11 

80.548 

446.99 

15.66 

414.66 

14.51 

417.94 

14.64 

89.981 

413.00 

13.54 

382.50 

12.49 

377.78 

12.37 

100.278 

468.11 

16.16 

469.20 

14.47 

449.23 

13.59 

111.715 

484.68 

25.36 

486.70 

24.74 

470.53 

122.090 

484.87 

25.49 

496.99 

24.24 

466.13 

13.93 

513.32 

13.68 

478.41 

12.94 

140.716 

532.16 

14.53 

538.50 

14.31 

512.85 

13.50 

499.67 

13.47 

149.761 

560.22 

14.26 

527.99 

13.43 

517.27 

13.14 

158.110 

563.44 

17.91 

550.33 

14.32 

531.06 

13.72 

527.40 

13.60 

169.488 

576.83 

15.59 

571.12 

14.29 

552.32 

180.974 

610.85 

14.90 

591.22 

14.30 

572.84 

13.88 

Table  39:  Evaporative  heat  transfer  coefficient  (W/m^K)  for  thermocouple  5,  C?  = 
1.5  and  Or  =  10.0-g. 


Bottom 

Outboard 

Top 

Inboard 

Qixi 

/ig 

mjm 

/ig 

CBII 

mSSMi 

h/Q 

rnsmi 

390.35 

675.21 

164.90 

859.31 

211.69 

1020.30 

257.08 

187.16 

847.80 

107.02 

949.04 

121.45 

30.875 

1301.98 

110.09 

646.08 

51.20 

804.54 

64.16 

880.23 

71.07 

41.022 

1433.55 

217.38 

681.79 

95.15 

810.19 

184.75 

881.20 

212.59 

50.181 

1262.11 

156.87 

794.65 

149.84 

852.69 

169.34 

61.273 

1268.49 

736.49 

92.79 

775.69 

100.82 

70.586 

1185.01 

82.94 

768.87 

89.55 

80.548 

1174.06 

709.13 

25.13 

736.56 

26.24 

89.981 

1144.59 

39.25 

686.98 

22.69 

705.12 

23.45 

100.278 

1139.51 

18.67 

709.81 

21.65 

728.73 

22.39 

111.715 

1166.11 

70.17 

32.11 

752.17 

32.76 

122.090 

1173.63 

66.98 

731.70 

30.01 

745.27 

30.65 

129.122 

1160.12 

622.32 

16.74 

736.16 

19.81 

747.97 

20.21 

140.716 

1181.13 

31.67 

636.29 

16.63 

754.5'r 

19.82 

763.99 

20.18 

149.761 

1174.68 

756.12 

19.19 

766.27 

19.52 

158.110 

1142.55 

29.90 

19.10 

750.00 

19.27 

169.488 

1127.02 

28.70 

615.08 

15.35 

739.08 

18.51 

747.10 

18.78 

180.974 

1106.86 

27.28 

619.64 

15.01 

752.48 

18.89 

755.01 

19.36 
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Table  40:  Thermal  resistance,  TVansported  heat  and  Input  heat  for  G  =  0.5  and 

=  0.01-g. 


iJo.  (K/W) 

Ai?th  (K/W) 

Qt  (W) 

AQ,  (W) 

Oto  (W) 

AQin  (W) 

2.4929 

1.30160 

5.829 

3.042 

11.254 

2.220 

3.3485 

1.29290 

7.791 

3.007 

15.494 

2.615 

4.8585 

1.80183 

8.360 

3.100 

20.777 

3.039 

4.5051 

1.25605 

11.528 

3.214 

25.045 

3.331 

Table  41:  Thermal  resistance,  Transported  heat  and  Input  heat  for  G  =  1.0  and 

=  0.01-g. 


i2th  (K/W) 

Ai2th  (K/W) 

Qt  (W) 

AQt  (W) 

Oin  (W) 

AQin  (W) 

0.4635 

0.21369 

7.091 

3.208 

10.903 

2.180 

0.3256 

0.06563 

16.948 

3.321 

mmifim 

3.053 

0.5154 

0.06955 

25.104 

3.338 

30.494 

3.666 

0.8716 

0.09011 

32.608 

3.358 

40.053 

4.192 

1.2130 

0.10852 

38.928 

3.475 

49.903 

4.707 

Table  42:  Thermal  resistance,  TVansported  heat  and  Input  heat  for  G  =  1.5  and 

=  0.01-g. 


Rtb  (K/W) 

Ai?th  (K/W) 

Qt  (W) 

AQ,  (W) 

Qm  (W) 

AQin  (W) 

0.4663 

0.16527 

8.666 

3.019 

10.124 

2.096 

0.2707 

0.05004 

17.406 

3.067 

20.024 

2.955 

0.1908 

0.02399 

27.874 

3.254 

29.462 

3.592 

0.1835 

0.01844 

36.886 

3.352 

39.808 

4.177 

0.1768 

0.01534 

44.497 

3.365 

50.142 

4.697 

0.1652 

0.01357 

54.705 

3.213 

60.577 

5.165 

0.1858 

0.01073 

63.351 

3.256 

69.754 

5.544 

0.2151 

0.01312 

72.570 

3.114 

80.619 

5.994 

0.2707 

82.626 

3.720 

89.645 

6.273 

0.3419 

0.01330 

92.707 

3.487 

101.692 

0.4430 

0.01617 

98.670 

3.544 

110.802 
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Table  43:  Thermal  resistance,  Transported  heat  and  Input  heat  for  G  =  0.5  and  Or 

=  10.0-g. 


fith  (K/W) 

AiJ,h  (K/W) 

Qt(W) 

AO,  (W) 

Oin  (W) 

(W) 

0.2555 

0.06309 

7.131 

1.677 

10.735 

2.174 

0.4580 

0.06988 

10.823 

1.621 

15.694 

2.625 

0.8389 

0.10052 

13.599 

1.621 

20.869 

3.041 

1.2670 

0.13577 

15.065 

1.611 

25.298 

3.361 

1.5288 

0.16352 

17.230 

1.841 

29.973 

3.661 

1.6983 

0.15601 

19.268 

1.768  ^ 

33.634 

3.899 

1.9696 

0.15523 

21.642 

1.704 

40.170 

4.220 

2.0622 

0.15873 

22.892 

1.761 

43.2851 

4.404 

Table  44:  Thermal  resistance,  Transported  heat  and  Input  heat  for  G  =  1.0  and 

=  10.0-g. 


i?th  (K/W) 

AB,h  (K/W) 

<3,{W) 

AQ,  (W) 

Oin  (W) 

AQto  (W) 

0.3529 

0.09036 

6.921 

1.721 

10.966 

2.196 

0.2324 

0.03002 

14.516 

1.762 

20.617 

3.007 

0.01719 

22.005 

1.800 

29.937 

0.1737 

0.01143 

29.894 

1.789 

40.167 

4.210 

0.1578 

0.00851 

37.977 

1.842 

0.1405 

0.00662 

60.373 

5.173 

0.1183 

0.00516 

55.452 

2.047 

70.243 

5.595 

0.0853 

0.00360 

62.975 

2.119 

80.156 

5.972 

0.00322 

70.110 

2.108 

88.516 

6.298 

0.0765 

0.00299 

77.254  ^ 

2.227 

98.32^ 

6.636 

0.0590 

0.00271 

89.159 

2.405 

0.0577 

0.00242 

119.747 

7.323 

toms 

106.054 

2.594 

131.622 

7.673 

0.2620 

0.00670 

110.573 

2.701 

139.843 

7.929 
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Table  45:  Thermal  resistance,  TVansported  heat  and  Input  heat  for  G  =  1.5  and  Or 

=  10.0-g. 


i?th  (K/W) 

Aiith  (K/W) 

0.  (W) 

AQ,  (W) 

Oto(W) 

AQto  (W) 

0.5400 

0.12865 

7.502 

1.770 

10.281 

2.107 

0.3435 

0.04257 

14.906 

1.797 

20.038 

2.951 

0.2629 

0.02075 

23.055 

1.742 

30.875 

3.668 

0.2292 

0.03527 

31.702 

4.285 

41.022 

4.256  . 

0.2130 

0.02785 

37.905 

4.301 

50.181 

4.705 

0.1696 

0.01528 

47.508 

3.429 

61.273 

5.203 

0.1493 

0.01258 

54.078 

3.460 

70.586 

5.601 

0.1282 

0.00479 

61.757 

2.043 

80.548 

5.972 

0.1304 

0.00454 

69.557 

2.133 

89.981 

6.312 

0.1178 

0.00386 

77.290 

2.168 

100.278 

0.1028 

0.00485 

88.473 

3.701 

111.715 

0.0992 

0.00446 

95.964 

3.769 

122.090 

0.0941 

0.00289 

101.478 

0.0854 

0.00262 

112.422 

2.654 

140.716 

0.0809 

0.00251 

119.096 

2.705 

149.761 

8.164 

0.0771 

0.00248 

124.457 

2.857 

158.110 

8.445 

0.0729 

0.00217 

131.473 

2.941 

169.488 

8.708 

0.0652 

0.00203 

140.110 

3.026 

180.974 

9.002 
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9.4  Steady-State  Temperature  Distributions 
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(Oo)j;  (Oo)j;  (Oo)j:  (Oo)j; 


(Oo)j:  (Oo)j:  (Oo)j;  (OoU 


in  =  10.0  w 

=  15.0  W  H-- 
in  =  20.0  W  -B- 
=  25.0  W  -x  - 


=  30.0  W  -A- 


=  2.0  g,  G  =  0.5:  (a) 


(Oo)j;  {Oo)x  iOo)x  {Oo)x 


n  =  10.0  w 

“  =  15.0  W  H- 
„  =  20.0W  -0-- 
“  =  25.0  W  -X-- 
“  =  30.0  W  -A- 
n  =  35.0  W  -X- 
"  =  40.0  W  -o-- 
“  =  45.0  W  H-  - 


:  10.0  g,G  =  0.5: 


(a) 


(Oo)^//  (Oo)j[/  (Oo)iX  (0o)j[/ 


(Oo)j:  (Oo)j:  (Oo)j:  iOoU 


iOoU  iOo)x  iOo)i  iOoU 


r(°C)  T(“C)  T(°C)  T(°C) 


(OoU  (oJx  (o.)x  (OoU 


=  10.0  w 

=  20.0  W  -H- 
=  30.0  W 
=  40.0  W  -x - 
=  50.0  W 
=  60.0  W 
=  70.0  W  -O-- 
=  80.0  W  - 
=  90.0  W  HE3- 
=  100.0  W  -X- 
=  110.0  W  -A- 
=  120.0  W  -)K- 
=  130.0  W  O-- 
=  140.0  W  -+-■ 
=  150.0  W 
=  160.0  W  -X- 
=  170.0  W  -A  - 


7.  =  8.0  g,  G  —  1.5: 


(a) 


(Oo)j;  (Oo)j;  {Oo)x  (OoU 


